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Abstract 

lO ' Multidimensional noncommutative Laplace transforms over octo- 

nions are studied. Theorems about direct and inverse transforms and 
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other properties of the Laplace transforms over the Cayley-Dickson 
I algebras are proved. Applications to partial differential equations 

r~| . including that of elliptic, parabolic and hyperbolic type are investi- 

. gated. Moreover, partial differential equations of higher order with 

real and complex coefficients and with variable coefficients with or 
without boundary conditions are considered. 

> 

^ : 1 Introduction. 

^ I The Laplace transform over the complex field is already classical and plays 

O ! very important role in mathematics including complex analysis and differen- 

tial equations [29l [121 123] . The classical Laplace transform is used frequently 
^ ■ for ordinary differential equations and also for partial differential equations 

^ I sufficiently simple to be resolved, for example, of two variables. But it meets 

substantial difficulties or does not work for general partial differential equa- 
tions even with constant coefficients especially for that of hyperbolic type. 

To overcome these drawbacks of the classical Laplace transform in the 
present paper more general noncommutative multiparameter transforms over 
Cayley-Dickson algebras are investigated. In the preceding paper a noncom- 
mutative analog of the classical Laplace transform over the Cayley-Dickson 
algebras was defined and investigated [H] . This paper is devoted to its gener- 
alizations for several real parameters and also variables in the Cayley-Dickson 
algebras. For this the preceding results of the author on holomorphic, that 
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is (super) differentiable functions, and meromorphic functions of the Cayley- 
Dickson numbers are used [T71[IS]. The super-differentiabihty of functions of 
Cayley-Dickson variables is stronger than the Frechet's differentiabihty. In 
those works also a noncommutative line integration was investigated. 

We remind that quaternions and operations over them had been first de- 
fined and investigated by W.R. Hamilton in 1843 0. Several years later 
on Cayley and Dickson had introduced generalizations of quaternions known 
now as the Cayley-Dickson algebras [21 El HH EH] . These algebras, especially 
quaternions and octonions, have found applications in physics. They were 
used by Maxwell, Yang and Mills while derivation of their equations, which 
they then have rewritten in the real form because of the insufficient devel- 
opment of mathematical analysis over such algebras in their time [H [71 [13] . 
This is important, because noncommutative gauge fields are widely used in 
theoretical physics [27] . 

Each Cayley-Dickson algebra Ar over the real field R has 2^" generators 
{io, ^1, ^2' -i} such that iq = 1, = —1 for each j = 1,2, ...,2^' — 1, 
'ij'ik = —ikij for every l</i;7^j<2^ — 1, where r > 1. The algebra ^r+i is 
formed from the preceding algebra At with the help of the so-called doubling 
procedure by generator i2r. In particular, Ai = C coincides with the field of 
complex numbers, ^2 = H is the skew field of quaternions, ^3 is the algebra 
of octonions, A^ is the algebra of sedenions. This means that a sequence of 
embeddings ... ^ Ar ^ ^r+i ^ ... exists. 

Generators of the Cayley-Dickson algebras have a natural physical mean- 
ing as generating operators of fermions. The skew field of quaternions is 
associative, and the algebra of octonions is alternative. The Cayley-Dickson 
algebra Ar is power associative, that is, z^^^ = z^z^ for each n, m G N 
and z ^ Ar. It is non-associative and non- alternative for each r > 4. A 
conjugation z* = z oi Cayley-Dickson numbers z G Ar is associated with the 
norm \z\^ = zz* = z*z. The octonion algebra has the multiplicative norm 
and is the division algebra. Cayley-Dickson algebras Ar with r > 4 are not 
division algebras and have not multiplicative norms. The conjugate of any 
Cayley-Dickson number z is given by the formula: 

(Ml) Z* ■.= ^* -7]\. 

The multiplication in Ar+i is defined by the following equation: 

(M2) (e + r/I)(7 + SI) = (^7 - h) + (5^ + ^7)1 
for each ^, r], 'j, S G Ar, z := ^ + r]l e Ar+i, C '■= 1 + ^ A+i- 

At the beginning of this article a multiparameter noncommutative trans- 
form is defined. Then new types of the direct and inverse noncommutative 
multiparameter transforms over the general Cayley-Dickson algebras are in- 
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vestigated, particularly, also over the quaternion skew field and the algebra 
of octonions. The transforms are considered in Ar spherical and Ar Carte- 
sian coordinates. At the same time specific features of the noncommutative 
multiparameter transforms are elucidated, for example, related with the fact 
that in the Cayley-Dickson algebra Ar there are 2*^ — 1 imaginary generators 
{ii, ...,'i2'--i} apart from one in the field of complex numbers such that the 
imaginary space in Ar has the dimension 2^ — 1. Theorems about properties 
of images and originals in conjunction with the operations of linear com- 
binations, differentiation, integration, shift and homothety are proved. An 
extension of the noncommutative multiparameter transforms for generalized 
functions is given. Formulas for noncommutative transforms of products and 
convolutions of functions are deduced. 

Thus this solves the problem of non-commutative mathematical analysis 
to develop the multiparameter Laplace transform over the Cayley-Dickson 
algebras. Moreover, an application of the noncommutative integral trans- 
forms for solutions of partial differential equations is described. It can serve 
as an effective means (tool) to solve partial differential equations with real 
or complex coefficients with or without boundary conditions and their sys- 
tems of different types. An algorithm is described which permits to write 
fundamental solutions and functions of Green's type. A moving boundary 
problem and partial differential equations with discontinuous coefficients are 
also studied with the use of the noncommutative transform. 

Moreover, a decomposition theorem of linear partial differential operators 
over the Cayley-Dickson algebras is proved. A relation between fundamental 
solutions of an initial and component operators is demonstrated. In con- 
junction with a line integration over the Cayley-Dickson algebras and the 
decomposition theorem of partial differential operators it permits to solve 
partial differential equations linear with constant and variable coefficients 
and non-linear as well as boundary problems (see also [19]). Certainly, this 
approach effectively encompasses systems of partial differential equations, 
because each function / with values in the Cayley-Dickson algebra is the 
sum of functions fjij, where each function fj is real- valued. 

All results of this paper are obtained for the first time. 

2 Multidimensional noncommutative integral 
transforms. 

1. Definitions. Transforms in Ar Cartesian coordinates. 
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Denote by Ar the Cayley-Dickson algebra, < r, which may be, in 
particular, H = ^2 the quaternion skew field or O = ^3 the octonion algebra. 
For unification of the notation we put ^ = R, ^1 = C. A function / : R" — > 
Ar we call a function-original, where 2 < r, n e N, if it fulfills the following 
conditions (1 — 5). 

(1) . The function f{t) is almost everywhere continuous on R" relative to 
the Lebesgue measure A„ on R". 

(2) . On each finite interval in R each function gj{tj) — f{ti, ...,tn) by tj 
with marked all other variables may have only a finite number of points of 
discontinuity of the first kind, where t — {ti, i„) e R", tj j — 1, n. 
Recall that a point i^o £ R is called a point of discontinuity of the first type, 
if there exist finite left and right limits lim„_>.„(,^u<ug g{u) =: g{uo - 0) e A 
and lim„^„o,„>„o g{u) -. g{uo + 0) e A^. 

(3) . Every partial function gj{tj) — f{ti, ...,tn) satisfies the Holder con- 
dition: \gj{tj -\- hj) — gj{tj)\ < Aj\hj\°'^ for each \hj\ < 5, where < a ^ < 1, 
Aj — const > 0, 5j > are constants for a given t — (ii,...,i„) e R", 
j — 1, n, everywhere on R" may be besides points of discontinuity of the 
first type. 

(4) . The function f{t) increases not faster, than the exponential function, 
that is there exist constants — const > v — (vi, i'„), a_i,ai e R, 
where Vj e {—1, 1} for every j — 1, n, such that 

\f{t) I < Cy exp((g^, t)) for each t e R" with tjVj > for each j — 1, n, 
qy = (t^ia^i, ...,t'„a^„); where 

(5) {x,y) :— YJj=iXjyj denotes the standard scalar product in R". 
Certainly for a bounded original / it is possible to take a_i = oi = 0. 
Each Cayley-Dickson number p e A we write in the form 

(6) p — Yl'f=Q Pjij-i where {io, ^2'--i} is the standard basis of gen- 
erators of Ar so that iQ — 1, Vj — —1 and igij — ij — ijio for each j > 0, 
ijik — —ikij for each j > and A; > with k ^ j, pj e R for each j. 

If there exists an integral 

(7) F"(p) - F-{p; C) := /r„ /(i)e-<^''*)-^di, 

then F"'{p) is called the noncommutative multiparameter (Laplace) trans- 
form at a point p e A of the function-original /(i), where C ~ Co = Ci^i + 
... + C2'--i^2'--i £ Ar is the parameter of an initial phase, Q e R for each 
j = 0, 1, 2^ - 1, C e A, ^ = 2^ - 1, dt = K{dt), 

(8) <p,t) ^ po{ti + ... + t2r-i) + Y^jJi Pjtjij, we also put 
(8.1) u{p,t;0 ^<p,t) + C. 

For vectors v,w E R" we shall consider a partial ordering 

(9) V ~< w if and only if Vj < Wj for each j — 1, ...,n and a k exists so 
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that Vk < Wk, 1 < k < n. 

2. Transforms in Ar spherical coordinates. 

Now we consider also the non-hnear function u — u(p,t;Q taking into 
account non commutativity of the Cayley-Dickson algebra Ar- Put 

(1) u{p, t) := u{p, t] C) := PoSi + M{p, t) + Co, where 

(2) M(p, t) = M(p, t; C) = (pisi + Ci)[ii cos(p2S2 + C2) + «2 sin(p2S2 + C2) 

C0S(P3S3 + Ca) + ■■■ + ^2'--2 Sin(p2S2 + C2) ■■■ Sin(p2'--2S2'--2 + C2'--2) COs(p2'--lS2'--l 

+C2'--i) + ^2'--i sin(p2S2 + C2)--- sin(p2'--2S2'--2 + C2'--2) sin(p2'--iS2'--i + C2'--l)] 
for the general Cayley-Dickson algebra with 2 < r < 00, 

(2.1) Sj :— Sj{n;t) :— tj + ... + tn for each j — 1, n = 2^ — 1, so 
that si — ti + ... + tn, Sn — tn. Morc generally, let 

(3) u{p,t) — u(p,t-X) — poSi + w{p,t) + Co, where w(p,t) is a locally 
analytic function, Re{w{p,t)) — for each p & Ar and t e R^"^"^, Re{z) :— 
{z -\- z)/2, z — z* denotes the conjugated number for z & Ar. Then the 
more general non-commutative multiparameter transform over Ar is defined 
by the formula: 

(4) F-{p- C) := V f{t) exp(-i.(p, C))dt 

for each Cayley-Dickson numbers p & Ar whenever this integral exists as the 
principal value of either Riemann or Lebesgue integral, n — 2^ — 1. This 
non-commutative multiparameter transform is in Ar spherical coordinates, 
when u(p, t; C) is given by Formulas (1, 2). 

At the same time the components pj of the number p and Q for C in 
u(p, t; C) we write in the p- and C-representations respectively such that 

(5) hj = {-hij + ij{2'- - 2)-^{-h + T,ZJi^ik{hil)})/2 for each j = 
1,2, ...,2^-1, 

(6) ho^{h + (2^ - 2)-'{-h + ElJi ik{hil)})/2, 

where 2 < r & N, h — hoio + ... + /i2'--i^2'--i £ Ar, hj e R for each j, 
il — ik — —ik for each k > 0, io — 1, h e Ar. Denote F^{p; () in more details 
hyJ'^if,u;p;C). 

Henceforth, the functions u{p, t; C) given by 1(8, 8.1) or (1, 2, 2.1) are used, 
if another form (3) is not specified. If for u{p, t; C) concrete formulas are not 
mentioned, it will be undermined, that the function u{p,t;C,) is given in Ar 
spherical coordinates by Expressions (1, 2, 2.1). If in Formulas 1(7) or (4) the 
integral is not by all, but only by tj{i), variables, where 1 < k < n, 

1 ^ i(l) < ■■■ < j{k) < n, then we denote a noncommutative transform by 
or jrfc;*^(i)'-'*^w(/,ii;p;C). If = l,...,j{k) = k, then we 
denote it shortly by F^{p;() or T^{f,u;p;C,). Henceforth, we take Cm = 
and tm — ^ and Pm — ^ for each 1 < m ^ 0(1), ■■■,i(^)} if something other 
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is not specified. 

3. Remark. The spherical Ar coordinates appear naturally from the 
following consideration of iterated exponents: 

(1) exp(H(piSi+Ci) exp(-i3(p2S2+C2) ^wi-iiipsSs+Cs))) = exp(ii(piSi+ 

Cl) exp(-(p2S2 + C2)(«3 C0S(P3S3 + (s) " «2 sin(p3S3 + (3)))) 

= exp(ii(pisi + Ci)(cos(p2S2 + C2) - sin(p2S2 + C2)(«3 cos(p3S3 + Cs) - 
i2 sin(p3S3 + C3)))) 

= exp((pisi + Ci)(«iCOs(p2S2 + C2) + i2sin(p2S2 + C2)cos(p3S3 + C3) + 
^3 sin(p2S2 + C2) sin(p3S3 + Ca)))- 

Consider the generator of the doubling procedure of the Cayley-Dickson 
algebra Ar+i from the Cayley-Dickson algebra Ar, such that iji2r — h^+j for 
each j — 0, 2^ — 1. We denote now the function M(p, t; () from Definition 
2 over Ar in more details by rM. 

Then by induction we write: 

(2) exp( r+iM{p,t;()) = exp{ rM{{iipi+,..+i2r-ip2r-i), {h, ...,t2r-2, {t2r-i+S2r)); 
(iiCi + ... + i2'--iC2'--i) exp(-i2r-+i(p2'-S2'- + C2O 

exp(- rM{{iiP2r+i + ...+i2r-lP2r+l-l), it2^+l, ■-, V+l-l); (^lC2'-+l + --+«2'--lC2'-+l-l)))}, 

where t = (ti, ...,tn), n = 'n{r + 1) = 2^+^ — 1, sj = Sj{n{r + for each 
i = 1, n(r + 1), since s^(n(r + 1); t) = + ... + tn(r+i) = Sm{n{r);t) + 
S2r{n{r + 1); t) for each m = 1, 2^ — 1. 

An image function can be written in the form 

(3)F„"(p;C):=E?:o'^,F„"^feC), 
where a function / is decomposed in the form 

(3.1) fit) = Y.U'i,f,{t), 

fj : R" — 7> R for each j = 0, 1, 2'' — 1, F^j(p; () denotes the image of the 
function-original fj. 

If an automorphism of the Cayley-Dickson algebra Ar is taken and instead 
of the standard generators {io, ^2'--i} new generators {A^"o, .... A^2'--i} ai'c 
used, this provides also M{p, t; Q = Mi\f(p, t; Q relative to new basic genera- 
tors, where 2 < r G N. In this more general case we denote by ^F^^p; () an 
image for an original f{t), or in more details we denote it by nJ^^H', u; p; (). 

Formulas 1(7) and 2(4) define the right multiparameter transform. Sym- 
metrically is defined a left multiparameter transform. They are related by 
conjugation and up to a sign of basic generators. For real valued originals 
they certainly coincide. Henceforward, only the right multiparameter trans- 
form is investigated. 

Particularly, if p = (po,Pi, 0, 0) and t = (ti, 0, 0), then the multi- 
parameter non-commutative Laplace transforms 1(7) and 2(4) reduce to the 



6 



complex case, with parameters Thus, the given above definitions over 

quaternions, octonions and general Cayley-Dickson algebras are justified. 
4. Theorem. // an original f{t) satisfies Conditions i(l — 4) and ai < 

a_i, then its image J-""- {f , u; p; () is Ar-holomorphic (that is locally analytic) 
by p in the domain {z G Ar '■ cti < Re{z) < ct-i}, as well as by ( ^ At, where 
1 < r G N, 2^^^ <2'^ — 1, the function u{p, t; () is given by i(8, 8.1) or 
^(1,2,2.1). 

Proof. At first consider the characteristic functions Xu^{t), where Xuit) = 
1 for each t G U, while xu{t) = for every t G R" \ t/, f/^ := {t G R" : 
Vjtj > Vj = 1, ...,«} is the domain in the Euclidean space R" for any v 
from §1. Therefore, 

(1) KiP;0 ■■= E[^=(,.l,...,^,„):^,l,...,^„e{-l,l}] lu. f {t) exp{-u{p,t; C))dt, 
since XniUyHUw) = for each v ^ w. Each integral /^^ f{t) exp(— t; C))dt 
is absolutely convergent for each p E Ar with the real part Oi < Re{p) < a_i, 
since it is majorized by the converging integral 

(2) \ Ju, f{t) exp{-u{p,t; C))dt\ < ... C^exp{-vi{w - ay^)yi - ... - 
Vn{w - ttyjyn - (o}dyi...dyn = C^e'^" n"=i Vj{w - a^J"\ 

where w = Re{p), since |e^| = exp{Re{z)) for each z E Arin view of Corollary 
3.3 [15]. While an integral, produced from the integral (1) differentiating by 
p converges also uniformly: 

(3) I / f{t)[dexp{-u{p,t;0)/dp].hdt\ 

< ... Cy\{ho{viyi + ...+Vnyn),hi{Viyi + ...+Vnyn),---,hn-^liVn-iyn-l+Vnyn),KVnyn)\ 

Jo Jo 

exp{-i;i(w - a^Jyi - ... - Vn{w - a^Jyn - Co}dyi...dyn 

n 

< |/i|ae-^«n(^ 

for each h E Ar, since each z E Ar can be written in the form z = \z \ exp(M), 
where \z\'^ = zz E [0, oo) C R, M G A, Re{M) := (M + M)/2 = in 
accordance with Proposition 3.2 [16j. In view of Equations 2(5,6): 

(4) fit) exp(-n(p, t; C))dt)/dp = and 

(5) dU^r. fit) expi-uip, t- C))dt)/dC = 0, while 

(6) \fuJit)[d^M-uip,t-C))/dC\.hdt\ < \h\J,^...J,^C,exp{-v^iw ~ 
a^Jyi - ... - Vniw - a^Jyn - Co]dyi...dyn = |/i|C^e"'^o n"=i Vjiw - a^^)~'^ 
for each h E Ar. In view of convergence of integrals given above (1 — 6) the 
multiparameter non-commutative transform -F"(p; C) is (super) differentiable 
by p and (, moreover, dF^ip;()/dp = and dF^ip;()/d( = in the con- 
sidered ip, C)-representation. In accordance with [T71 [TB] a function gip) 
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is locally analytic by p in an open domain U in the Cayley-Dickson alge- 
bra At, 2 < r, if and only if it is (super) differentiable by p, in another 
words ^r-holomorphic. Thus, F^{p;() is ^^-holoniorphic hj p E Ar with 
ai < Re{p) < a_i and C £ Ar due to Theorem 2.6 [T8] . 

4.1. Corollary. Lei suppositions of Theorem 4 be satisfied. Then the 
image J^"'{f,u;p;() with u = u{p,t;() given by 2{1,2) has the following pe- 
riodicity properties: 

(1) C + l^ij) = J^^if, u;p; Q for each j = 1, n and (5 G 27rZ; 

(2) F^{f,u-p^-C) = {-lYF^{f,u-p'-e) for each] = so that 
Co = Co o.nd Cj = -Cj, Cj+i = TT + Cf+i, Ci = C.^ for each j ands^ j + l, 
while either pj = —p^ and p] = pf for each I ^ j with k = 2 or p^ = p^ and 
f{t) is an even function with n = 2 by the Sj = {tj + ... + tn) variable or an 
odd function by Sj = {tj + ... + t^) with n = 1; 

(3) J^^if, u-p;C + vr^i) = u; p; () . 

Proof. In accordance with Theorem 4 the image J^'"{f, u]p; C) exists for 
each p G Wj := {z E Ar : ai < Re{z) < a_i} and ( G Ar, where 1 < r. 
Then from the 27r periodicity of sine and cosine functions the first statement 
follows. From sin(— 0) = — sin(0), cos(0) = cos(— 0), sin(7r + 0) = — sin(0), 
cos(0 + 7r) = — cos(0) we get that cos{pjSj + (j) = cos{—pjSj + (j), sm{pjSj + 
Q) cos{pj+iSj+i + C}+i) = {-sm{-pjSj + C|))(-cos(pj+iSj+i + C|+i)) and 
sm{pjSj+Cj) sin(pj+iSj+i + C]+i) = (- sm{-pjSj + Cf)){- sin(pj+iSj+i+C|+i)) 
for each t G R". On the other hand, either p^ = —p'j and pj = pf for each 
/ 7^ j > 1 with K = 2 orp^ = p2 and /(ti, Sj_i+Sj, -Sj-Sj+i,tj+i, = 
(— Sj^i — Sj, Sj — Sj+i, tj+i, ...,t„) is an even with k = 2 or odd with 
K, = 1 function by the sj = {tj + ... variable for each t = (ti, ...,tn) G R", 
where tj = Sj — Sj^i for j = l,...,n, Sn+i = Sn+i{n;t) = 0. From this 
and Formulas 2(1,2,4) the second and the third statements of this corollary 
follow. 

5. Remark. For a subset U in Ar we put 7Cs.p.t{U) := {u : z E U, z = 
EvGb^^vV, u = WsS + Wpp} for each s 7^ p G b, where t := Eveb\{s,p} ^^'vV G 
^r,s,p := {z e Ar : z = I^vGb^vV, Ws = Wp = 0, Wy E H E b}, where 
b := {zo, ^1, '^2'--i} is the family of standard generators of the Cayley- 
Dickson algebra Ar. That is, geometrically 7Cs,p,t{U) means the projection 
on the complex plane Cj^p of the intersection U with the plane 7fs,p,t 3 t, 
Cs,p := {as + bp : a,b E R}, since sp* G 6 := b \ {1}. Recall that in §§2.5-7 
[16] for each continuous function f : U Ar it was defined the operator / 
by each variable z E Ar- For the non-commutative integral transformations 
consider, for example, the left algorithm of calculations of integrals. 

A Hausdorff topological space X is said to be n-connected for n > if 
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each continuous map f : ^ X from the /c-dimensional real unit sphere 
into X has a continuous extension over R'*^"''-'- for each k < n (see also |28]). 
A 1-connected space is also said to be simply connected. 

It is supposed further, that a domain U in Ar has the property that 
U is (2'" — l)-connected; iTs^p^tiU) is simply connected in C for each k = 
0, 1, 2'""-'^, s = i2k, P = i2k+i, t G ^r,s,p and u G Cs,p, for which there exists 
z = u + teU. 

6. Theorem. // a function f{t) is an original (see Definition 1), 
such that j\fF^ {p; () is its image multiparameter non- commutative transform, 
where the functions f and F" are written in the forms given by 3{3,3.1), 
/(R") C Ar over the Cayley-Dickson algebra Ar, where 1 < r G N, 2^^^ < 
n<T' -I. 

Then at each point t, where f{t) satisfies the Holder condition the equality 
is accomplished : 

/NnOO fN\00 
](-([(2vriV0-W ]NF:{a + p-C) 
-NnOO J-NlOO 

exp{w(a + p, t- 0})-)dp} =: NF^ia + p; C),u, t; C), 

where either u{p, t; () =< p,t) + ( or u{p, t; Q = poSi + M^ip, t; () + (q (see 
§§i and 2), the integrals are taken along the straight lines p{Tj) = NjTj G Ar, 
Tj G R for each j = l,...,n; ai < Re{p) = a < a_i and this integral is 
understood in the sense of the principal value, t = {ti,...,tn) G R", dp = 
{...{{d[p,N,])d[p2N2])...)d[pnN,,]. 

Proof. In Integral (1) an integrand 7]{p)dp certainly corresponds to the 
iterated integral as {...{ri{p)d[piNi])...)d[pnNn], where p = piNi + ... +pnNn, 
Pi, ...,Pn G R. Using Decomposition 3(3.1) of a function / it is sufficient to 
consider the inverse transformation of the real valued function fj, which we 
denote for simplicity by /. We put 

ivF.",(p;C):= / fjit)exp{-uip,t;C))dt. 

If ?7 is a holomorphic function of the Cayley-Dickson variable, then locally 
in a simply connected domain U in each ball B{Ar, Zq, R) with the center at 
zq of radius R > contained in the interior Int{U) of the domain U there is 
accomplished the equality 

id[!:^via + OdC]/dz).l = r^ia + z), 
where the integral depends only on an initial zo and a final z points of a rec- 
tifiable path in B{Ar, Zq, R), a E H (see also Theorem 2.14 [18j). Therefore, 
along the straight line iV^-R the restriction of the antiderivative has the form 
ll via + NjTj)dTj, since 
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(2) i:-=Je, via + OdC = Jl fj{a + N,T,).N,dT„ 

where dr]{a + z)/d9 = {dT]{a + z)/dz).Nj for the (super) differentiable by 
z & U function ri{z), when z = ONj, 6 E H.. For the chosen branch of the 
hne integral specified by the left algorithm this antiderivative is unique up 
to a constant from Ar with the given ^-representation u of the function rj 
fl6[ \T7\ [18]. On the other hand, for analytic functions with real expansion 
coefficients in their power series non- commutative integrals specified by left 
or right algorithms along straight lines coincide with usual Riemann integrals 
by the corresponding variables. The functions sin{z), cos(z) and partici- 
pating in the multiparameter non-commutative transform are analytic with 
real expansion coefficients in their series by powers of z G Ar- 

Using Formula 4(1) we reduce the consideration to XuS'^)fif) instead of 
/(t). By symmetry properties of such domains and integrals and utilizing 
change of variables it is sufficient to consider with f = (1, 1). In this 
case /j^„ for the direct multiparameter non-commutative transform 1(7) and 
2(4) reduces to /q°°.../q°°. Therefore, we consider in this proof below the 
domain only. Using Formulas 3(3,3.1) and 2(1,2,2.1) we mention 

that any real algebra with generators A^o = 1, ^fc and Nj with 1 < A; 7^ j 
is isomorphic with the quaternion skew field H, since Re{NjNk) = and 
\Nj\ = 1, \Nk\ = 1 and \NjNk\ = 1. Then exp(a + M/3)exp(7 + Mu) = 
exp((a+7)+M(/3+ci;)) for each real numbers a, (3, 7, 6 and a purely imaginary 
Cayley-Dickson number M. 

The octonion algebra O is alternative, while the real field R is the center 
of the Cayley-Dickson algebra Ar- We consider the integral 

(3) g,it) := [i2nNr,y' jXb]i-m^N,)-' iXb] NF^/a+p; () exp{uia+ 
p,t-X)})-)dp 

for each positive value of the parameter < 6 < 00. With the help of 
generators of the Cayley-Dickson algebra Ar and the Fubini Theorem for 
real valued components of the function the integral can be written in the 
form: 

(4) (7fe(t) = [(27riV„)-W rfr„/ ](...([(27riVi)-W dn ] 

Jo J-N„b Jo J-Nib 

f{T) exp{-UAr(a + p, t; C)} exp{MAr(a + p, r; C)})--)dp, 

since the integral J^_^ ^ /{t) exp{—UN{a + p,T; ()}dT for any marked < 
6 < (a_i — ai)/3 is uniformly converging relative to p in the domain ai + 6 < 
Re{p) < a_i — 5 m Ar (see also Proposition 2.18 [T8]). 

If take marked tk for each k ^ j and S = Nj for some j > 1 in Lemma 2. 17 
|18j considering the variable tj, then with a suitable (R-linear) automorphism 
V of the Cayley-Dickson algebra Ar an expression for v(M(p, t; ()) simplifies 
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like in the complex case with Ck '■— R®Ri^ for a purely imaginary Cayley- 
Dickson number K, \K\ — 1, instead of C := R © Rii, where v{x) — x for 
each real number x e R. But each equality a = ^ in is equivalent to 
v(q;) = v(^). Then 

(5) Re[(NjN,)(NjNiY] = Re(N,N^) = 5,,i for each q,l. 

If = Eo<i<n;i^j (^i^h = Eo<i<n;i^j M with j > 1 and real num- 
bers ai,/3i e R for each I, then 

(6) Re[{NjS^){NjN^y] = Re[S^{N^y] = EiO^iPi- 

The latter identity can be applied to either 5''^ = Mk+i{pk+iNk+i + ... + 

PnNn, {tk+l, -, tn); Ck+lNk+1 + ... + Cn^n) and A^'^ = Mk+i{pk+lNk+i + ... + 
PnNn, {Tk+l, Tn); Cfe+1-/Vfe+1 + ■■■ + CnNn), Or S'' = {pk+ltk+1 + Ck+l)Nk+l + 
... + {Pntn + Cn)Nn and N'' = (pfe+iTfe+i + Cfe+l)-/Vfe+l + ■■■ + {PnT„ + Cn)Nn, 

where 

(7) Mk+i{pk+lNk+i + ... + PnNn, {tk+l, tn); Ck+l^k+l + ... + Cn^n) = 
iPk+lSl,k+l + Cfe+l)[-^fe+l COs(pfc+2S2,fe+l + Cfe+2) + ■■■ 

+A^„Sin(pfe+2S2,fe+l + Cfe+2)-- Sm{pnSn-k,k+l + Cn)], 

(8) Sj^k+1 = Sj^k+i{n;t) = tk+j + ...+tn = Sk+j{n;t) for each j = 1, 

Sn-k,k+l — Sn-k,k+l{'>T';t) — tn. 

We take the limit of gb{t) when b tends to the infinity. Evidently, Sk{n; r) — 
Sj{n;T) — Sk{j — 1;t) — Tk + ... + tj-i for each 1 < k < j < n. By our 
convention Sk{n; r) = Si{n; r) for k < 1, while Sk{n; t) — ior k > n. Put 

(9) Un,j{po+PjNj + ... +PnNn, {tj, Tn); Co + Cj^j + ... + Cn-^n) = Co + 
PoSij + Mj{pjNj + ... +PnNn, {tj, T„); Co + CjNj + ... + CnNn) 

for un given by 2(1, 2, 2.1), where Mj is prescribed by (7), Skj — Skj{n; r); 

(10) Un,j{po + PjNj + ... + PnNn, {tj, ...,r„);Co + CjNj + ... + CnNn) = 
Co + P0Si,j + T,k=j{Pk'Tk + Ck)Nk 

for u — Un given by 1(8, 8.1). For j > 1 the parameter Co for u — un given 
by 1(8, 8.1) or 2(1, 2, 2.1) can be taken equal to zero. 

When ii, f^+i, tn andpi, variables are marked, 

we take the parameter 

C^' := CKVjNj + - + PnNn, {Tj, Tn); Co + CjNj + ... + CnNn) := (Co + 

CjNj + ... + CnNn) + (a + po)sj+i +Pj+iSj+iNj+i + ...+pnSnNn for u{p, t; C) 
given by Formulas 2(1, 2, 2.1) or 

C^' := C'{PjNj + ... + PnNn), {Tj, Tn); Co + CjNj + ... + CnNn) := (Co + 
CjNj + ... + CnNn) + {a+Po)Sj+i + Pj+iTj+iNj+i + ... +PnTnNn for u{p, T; C) 

described in 1(8,8.1). Then the integral operator 

\imb^^[{27r Nj)~^ /q°° drj j!^^.^]...{dpjNj) (see also Formula (4) above) applied 

to the function f{ti, tj^i, Tj, r„) exp{-UN,j{a+po+PjNj+...+PnNn, {tj, tn); Co+ 
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... + CnNn)} with the parameter instead of ( treated by Theorems 2.19 
and 3.15 [18] gives the inversion formula corresponding to the real variable 
tj for f{t) and to the Cayley-Dickson variable PqNq +PjNj restricted on the 
complex plane Cjy. = RSRA^^-, since d{Tj + c) = rfr,- for each (real) constant 
c. After integrations with j = 1, k with the help of Formulas (6 — 10) and 
3(1,2) we get the following: 

(11) lira g,{t) = Re[{2nN„)-' drj ]{...{[{2nN,+,)-' dn+i ] 

b^co Jo J-Nnoo Jo J-Nk+lOO 

/(ti, tfc, Tfc+i, r„) exp{-UN,k+i{{a+Po+Pk+iNk+i+...+PnNn), (4+i, t„); 

(Co + Cfc+l^fc+1 + ••• + CnNn))} exp{uN,k+l{{a + Po+ Pk+lNk+l + ... + PnNn), 

(rfc+1, ...,r„); (Co + Cfc+iA^fc+i + ••• + C,nNn))})-:)dp. 

Moreover, Re{fg) = fg for each q and in (11) the function f = fq stands for 
some marked q in accordance with Decompositions 3(3,3.1) and the begin- 
ning of this proof. 

Mention, that the algebra algji{Nj, Nk, Ni) over the real field with three 
generators Nj, and A*"; is alternative. The product Nj^Ni of two generators 
is also the corresponding generator (— 1)^*^'^''^^!^ with the definite number 
m = m{k,l) and the sign multiplier (— 1)^*^^''\ where ^{k,l) G {0,1}. On 
the other hand, Nkj^[Nj{Nj{Nk,^Ni))] = Nk^{Nk^Ni). We use decompositions 
(7 — 10) and take k2 = I due to Formula (11), where Re stands on the right 
side of the equality, since Re^N^Ni) = and Re[Nj{Nj{NkNi))] = for each 
k ^ I. Thus the repeated application of this procedure by j = 1, 2, n leads 
to Formula (1) of this theorem. 

6.1. Corollary. If the conditions of Theorem 6 are satisfied, then 

(1) f{t) = {27r)-- f F:{a + p;C)exp{u{a + p,t;C)}dpi...dpn 
= (^")-i(;vF:(a + p;C),u,t;C). 



Proof. Each algebra algii{Nj, Nj^, Ni) is alternative. Therefore, in accor- 
dance with §6 and Formulas 1(8, 8.1) and 2(1 — 4) for each non-commutative 
integral given by the left algorithm we get 

fNjb 

(2) J_^Jfir) exp{-UN{a+p,t; ()}] exp{M^(a +p, r; C)}d{pjNj) 

J2 Nj[Nj{ [Nifi{T)exp{-UN{a + p,t;C)}]exp{uN{a + p,T;C)}dp,)] 
1=0 ''-^^^ 
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= / [f{r)exp{-UN{a+p,t;C)}]ew{uN{a+p,T;C)}dpj 

J —h 

for each j = 1, since the real field is the center of the Cayley-Dickson 
algebra Ar-, while the functions sin and cos are analytic with real expansion 
coefficients. Thus 

(3) g,{t) = (27r)-"[/o°^rfr„/^,](...([/o-dTi/!,]/(r)exp{-«^(a+p,t;C)} 

exp{MAr(a + p, t; ()})... )dpi...dpn, 
hence taking the limit with b tending to the infinity implies, that the non- 
commutative iterated (multiple) integral in Formula 6(1) reduces to the prin- 
cipal value of the usual integral by real variables (ri, ...,Tn) and (pi, ...,Pn) 
6.1(1). 

7. Theorem. An original f{t) with /(R") C Ar over the Cayley- 
Dickson algebra Ar with 1 < r E 'N is completely defined by its image 

nF^{p; Q up to values at points of discontinuity, where the function u{p,t; Q 
is given by i(8,8.1) or 2(1,2,2.1). 

Proof. Due to Corollary 6.1 the value f(t) at each point t of continuity of 
f{t) has the expression throughout nF^{p; Q prescribed by Formula 6.1(1). 
Moreover, values of the original at points of discontinuity do not infiuence on 
the image nF^{p', C)j since on each bounded interval in R by each variable 
tj a number of points of discontinuity is finite and by our supposition above 
the original function f{t) is A„-almost everywhere on R" continuous. 

8. Theorem. Suppose that a function iqF^{p; Q is analytic by the 
variable p G Ar in a domain W := {p E A,- : ai < Re{p) < a_i}, where 
2 < r G N, T'-^ < n < 2'' - 1, /(R") C A, either u{p, t; () =< p,t) + C or 
u{p,t; () := PqSi + M{p,t; + (q (^see §§i and 2). Let nF^{P'X) be written 
in the form nF^{P',0 = JV-f"u + nF^'^{P',C), where jvi^"'°(p;C) 
holomorphic by p in the domain ai < Re{p). Let also jv-P"^''^(p; C) be holo- 
morphic by p in the domain Re{p) < a_i. Moreover, for each a > ai and 
b < a_i there exist constants Ca > 0, Cb > and > and > such that 

(1) I NF^'^ip; 01 < Caexp(-e„|p|) for each p E Ar with Re{p) > a, 

(2) I nF^'\p;C)\ < C6exp(-e6|p|) for each p G Ar with Re{p) < b, the 
integral 

(3) J!^nToo ■■■ I-nToo nF^^'^{w + p; C)dp converges absolutely for k = and 
k = 1 and each ai < w < a_i. Then nF^{w + p;C) ^■s the image of the 
function 

]{...{[{2n)-'N, / ] ^F:{w+p;C)exp{u{w+p,t;C)})...)dp 
^(J^)-\^F:(w+p;0,u,t;0. 
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Proof. For the function ^F^'^lp; () we consider the substitution of the 
variable p = —g, — a_i < Re{g). Thus the proof reduces to the consideration 

of nF:''{w + p;C). 

An integration by dp in the iterated integral (4) is treated as in §6. Take 
marked values of variables pi, ...,p„ and ti, tj+i, 

where = Sk{n; r) for each k = 1, ...,n (see §6 also). For a given parameter 

:= (Co + (jNj + ... + CnNn) + {W + Po)Sj+l + Pj+lSj^iNj+i + ... + PnSnNn 

for u{p,T-X) prescribed by Formulas 2(1,2,2.1) or := (^o + Cj^j + ••• + 
CnNn) + (w + Po)sj+i + pj+iTj^iNj^i + ... + PnTnNn for u{p,t] Q giveu by 
1(8,8.1) instead of C and any non-zero Cayley-Dickson number /3 G we 
have lim,^.^oo[/3r,- + C'']/[/3r, +(] = !• 

For any locally ^-analytic function g{z) in a domain U satisfying con- 
ditions of §5 the homotopy theorem for a non- commutative line integral 
over r, is satisfied (see [TTl [IS]). In particular if U contains the 

straight line w + 'RNj and the path 7j (tj) := C"* + tjNj, then g{z)dz = 

J^^ g{w + z)dz, when g{z) — )■ while \z\ tends to the infinity, since is a 
finite number (see Lemma 2.23 in [H]). We apply this to the integrand in 
Formula (4), since i\fF^{w + p; () is locally analytic by p in accordance with 
Theorem 4 and Conditions (1,2) are satisfied. 

Then the integral operator [(27rA'j)~^ I-n^oo\ ^^e j'-th step with the help 
of Theorems 2.22 and 3.16 [18] gives the inversion formula corresponding to 
the real parameter tj for f{t) and to the Cayley-Dickson variable poNo+pjNj 
which is restricted on the complex plane C^. = R © (see also Formulas 
6(4, 11) above). Therefore, an application of this procedure by j = 1, 2, n 
as in §6 implies Formula (4) of this theorem. 

Thus there exist originals and for functions ^F^'^^p; () and ^F^'^^p] () 
with a choice of w G R in the common domain Oi < Re{p) < a_i. Then 
/ = /° + is the original for nF^^p'X) due to the distributivity of the 
multiplication in the Cayley-Dickson algebra Ar leading to the additivity of 
the considered integral operator in Formula (4). 

8.1. Corollary. Let the conditions of Theorem 8 he satisfied, then 

(1) f{t) = (27r)-" / nF:{w +p-Oe^^{u{w + p,t-0}dpi...dpn 
= (.F")-i(^F:(«;+p;C),«,t;C). 

Proof. In accordance with §§6 and 6.1 each non- commutative integral 
given by the left algorithm reduces to the principal value of the usual integral 



14 



by the corresponding real variable: 

/Njoo 
nF:{w + p- C) exp{M(«; + p, t; 0}d{p,N,) 
-NjOo 

/oo 
nF::{w + p; C) exp{u{w + p, t; Qjdpj 
-oo 

for each j = 1, ...,n. Thus Formula 8(4) with the non-commutative iterated 
(multiple) integral reduces to Formula 8.1(1) with the principal value of the 
usual integral by real variables (pi, ...,p„). 

9. Note. In Theorem 8 Conditions (1, 2) can be replaced on 

(1) lim„^oo suppg(;^^^j ||^(p)|| = 0, 

where Ci?(„) := {z E Ar : \z\ = R{n),ai < Re{z) < a_i} is a sequence of 
intersections of spheres with a domain W, where R{n) < R{n + 1) for each 
n, lim„_>.oo R{n) = oo. Indeed, this condition leads to the accomplishment 
of the Ar analog of the Jordan Lemma for each r > 2 (see also Lemma 2.23 
and Remark 2.24 p^). 

Subsequent properties of quaternion, octonion and general Ar multipa- 
rameter non-commutative analogs of the Laplace transform are considered 
below. We denote by 

(2) Wf = {peAr: a,{f) < Reip) < a_i(/)} a domain of ^F^(p;C) 
by the p variable, where ai = ai(/) and a_i = a_i(/) are as in §1. For an 
original 

(3) f{t)xu,,...At) we put Wf = {peAr: ai{f) < Re{p)}, 

that is a_i = oo. Cases may be, when either the left hyperplane Re{p) = ai 
or the right hyperplane Re{p) = a_i is (or both are) included in Wj. It may 
also happen that a domain reduces to the hyperplane Wj = {p : Re{p) = 
ai = a^i}. 

10. Proposition. If images nF^{p'X) o.^'^d nG^^p'X) 

of functions- originals f{t) and g(t) exist in domains Wf and Wg with values 
in At, where the function u{p,t;Q is given by i(8,8.1) or ^(1,2,2.1), then 
for each a, P E Ar in the case A2 = H; as well as f and g with values in R 
and each a, (3 E Ar or f and g with values in Ar and each a, /3 G R m the 
case of Ar with r > 3; the function a nFu{p] C) + nGu{p', C) is the image 
of the function otf{t) + (ig{t) in a domain Wf fl Wg. 

Proof. Since the transforms nF^{p',C) and nG^^p'X) exist, then the 
integral 

/ {af{t) + Pg{t))exp{-u{p,t;C))dt= f af{t) exp{-u{p,t; C))dt 
+ [ (5g{t)exY>{-u{p,t-C))dt 
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converges in the domain 

^f^^g = {P e Ar : max(ai(/),ai(^)) < Re{p) < min(a_i(/), a„i(^))}. 
We have t G R*^, 2''^^ < n < 2''' — 1, while R is the center of the Cayley- 
Dickson algebra Ar- The quaternion skew field H is associative. Thus, under 
the imposed conditions the constants a, /3 can be carried out outside integrals. 

11. Theorem. Let a = const > 0, let also F"{p;() be an image of an 
original function f{t) with either u =< p, t)+( oru given by Formulas 2{1, 2) 
over the Cayley- Dicks on algebra Ar with 2 < r < oo, T'^^ < n <2^' — \. 
Then an image F^{p/a\ C)/*^" of the function f{at) exists. 

Proof. Since pjSj+(j = Pj{s'j/a)+(j = {pj/a)s'j+Q for each j = l,...,n, 
where Sja = s'j, Sj = Sj{n;t), s'j = Sj{n;T), tj = atj for each j = l,...,n. 
Then changing of these variables implies: 

/r" f{at)e-''^P'^'<^dt = /r„ /(r)e-"(P'^/"'^)dr/a" = 
due to the fact that the real filed R is the center Z{Ar) of the Cayley-Dickson 
algebra Ar- 

12. Theorem. Let f{t) be a function- original on the domain f/i,...,i such 
that df{t)/dtk also for k = j — 1 and k = j satisfies Conditions i(l — 4). 
Suppose that u{p,t-X) is given by ^(1,2,2.1) or i(8,8.1) over the Cayley- 
Dickson algebra Ar with 2 < r < oo, 2^^^ <n<2^ — \. Then 

(1) 7\{df{t)idt,)xij, At),u-p-o = -r^-'-'''{f{t)xu.,,,,At'),u{p,t^-,cy,p-,o 

+[po + jzvk%,]r''Uit)xu,_,it),u-p- 

fc=i 

in the Ar spherical coordinates or 

(1.1) 7-{{df{t)/dt,)xu._,{t),u-p-c) = -:F--'-''\f{t)xu._,{t^)Mv.t^-X):p-X) 
+bo + P,Se,]-^"(/(t)xi/,.....,(t),M;p;C) 

in the Ar Cartesian coordinates in a domain W = {p E Ar '■ max(ai(/), ai{df /dtj)) < 
Re{p)}, where P := (ti, ■■■,tn : tj = 0), Sg^ = —d/d(k for each k > 1. 
Proof. Certainly, 

(2) dfitis))/ds, = df{t)/dti and 

(2.1) df{t)/dt, = Tl=i{df{t{s))/ds,){ds,/dt,) = Ei=,df{t{s))/dsk 
for each j = 2, ...,n, since tj = Sj — Sj^i, ti = Si — S2, where Sj = Sj{n;t), 
Sn+i = for each Z > 1. From Formulas 30(6, 7) [18] we have the equality in 
the Ar spherical coordinates: 

(3) dexp{-u{p,t;C))/dsj = -poSij exp{-u{p,t; C))-pjSe^ exp{-u{p,t; 
since 

exp{-u{p, t; 0) = exp{-poSi - Co} exp{-M{p, t; C)), 
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(9exp(-j9oSi - Co)/dsj = -poSij exp(-poSi - Co), 

d[cos{pjSj + Q) - sm{pjSj + Cj)'''j]/9sj = dexp{-{pjSj + Cj)ij)/dsj 

= -pjijexp{-{pjSj + Cj)ij) = -pjexp{-{pjSj + Q - 7r/2)ij) 

= -vAco^ijPjSj + Ci - 7r/2) - sin(pjSj + Q - 7r/2)ij] = 

-Pj^e, [cos{pjSj + Q- sm{pjSj + Qij], 

since Sj and Sk are real independent variables for each k ^ j, where Sj^k = 
for j k, while 6jj = 1, 

(3.1) SeJcos(pjSj + Ci) - sin(pjSj + Cj)ij] = 
-d[cos{pjSj + Cj) - sm{pjSj + Cj)ij]/dCj 

= [cos{pjSj + Cj - 7r/2) - sin(pjSj + Cj - 7r/2)ij]. 

In the Cartesian coordinates we take tj instead of Sj in (3.1). If (j){z) 
is a differentiable function by Zj for each j, cj) : Ar ^ Ar, zj = pjtj + Q, then 

(3.2) 9exp(-</.(z))/9(gt,) = -q[dexp{0/d^]k=^Ud<Piz)/dz,)p, 

= —qpj{—dexp{—(f){z))/dQ) = — p^Sqe^ exp(— 0(z)), where either q = I 
or q = —1, since dzj/dQ = 1. 
That is 

(3.3) Sg_, exp(— Zfc((/)fc + Cfc)) = for each j ^ k > 1 and any positive 
number x > 0, 

(3.4) Sl^ exp{-ij{(f)j + Q)) = exp{-ij{(j)j + Q - x7t/2)) and 
S^e,. exp(-ij(0j + Cj)) = exp{-ij{^j + Cj + xn/2)) 

for each non-negative real number x >0, (pk and Cfc ^ R-, where Se^ = Sej (Cj), 
the zero power S^^ = / is the unit operator; 

(3.5) Sge,e-"(P'*'^) = e-P«"i-^« 

r/[io5j- 1 cos(piSi+Ci) + (l-5j,i)ij_i sin(piSi+Ci)... cos(pjSj+Cj)+{Efc=^-^ 4 sin(piSi + 
Ci)... cos{pk+iSk+i + Cfc+i)} + ^2--! sin(pisi + Ci)... sin(p2' -iS2'--i + C2'--i)] 
in the Ar spherical coordinates, where either g=lorg = — 1 and 

(3.6) 7;-e(0):=e(0-W2) 

for any function C.{Cj) and any real number x G R, where j > 1. Then in 
accordance with Formula (3.2) we have: 

(3.7) Sqe. exp(-M(p, t; C)) = 

= [En=i ElZlm^))'q^Mi^)r-'-'M]k=-u(pM) 
for u{p,t;C) given by Formulas 1(8,8.1) in the Ar Cartesian coordinates, 
where either g = lorg = — 1. 

The integration by parts theorem (Theorem 2 in §11.2.6 on p. 228 |lUj ) 
states: if a < 6 and two functions / and g are Riemann integrable on the 
segment [a, b], F{x) = A + fit)dt and G{x) = B + git)dt, where A and 
B are two real constants, then jj^ F{x)g{x)dx = F{x)G{x)\''^ — jj^ f{x)G{x)dx. 
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Therefore, the integration by parts gives 

(4) {df{t)/dtj)exp{-u{p,t;0)dt, = f{t) exp{-u{p,t-X))\t^Zo 



[f{t){dexp{-u{p,t;C))/dt,)]dt,. 

Using the change of variables t ^ s with the unit Jacobian d(ti, ...,tn)/d{si, s„) 
and applying the Fubini's theorem componentwise to fjij we infer: 

(5) / idfit)/dt,) exp(-w(p, t; C))dt = I {df{t)/dt,) exp(-w(p, t; C))ds 

poo poo poo 

= ... [ {df{t)/dtj)exp{-u{p,t;C))ds,]df 

Jo Jo Jsj^i 

poo poo 

= -[/ ... / f{tnexp{-u{p,t^;C))dt^] 
Jo Jo 

3 poo poo 

+ [Po + J2Pk^^k] - f{t)exp{-u{p,t;C))dt 
f— ^ Jo Jo 



k=l 

in the Ar spherical coordinates, or 

(5.1) / idfit)/dtj)expi-uip,t;C))dt 

poo poo 

= -[/ ... / f{tnexp{-u{p,t^;C))dt^] 
Jo Jo 

poo poo 

+\po+PjSe,] - f{t)exp{-u{p,t-X))dt 

in the Cartesian coordinates, since (?exp(—(poSi+Co))/f^^i = — Po exp(— (poSi+ 
Co)) for each I < j <n. This gives Formula (1), where 

(6) r'-''''{mxu,_„u{p,t^;0;P;0= ■■■ mexp{-uip,t^;0)dt' 

Jo Jo 

poo poo poo poo 

= / dti... / dtj^i / dtj+i... / {dtn)f{t^)exp{-u{p,t^;C)) 
Jo Jo Jo Jo 

is the non-commutative transform by P = (ti, 0, tj+i, ■■■,tn). 

12.1. Remark. Shift operators of the form C{x + 0) = cxp{(f)d/dx)^{x) 
in real variables arc also frequently used in the class of infinite diffcrcntiablc 
functions with converging Taylor series expansion in the corresponding do- 
main. 

It is possible to use also the following convention. One can put cos(0i + 
Ci) = cos(0i + Ci)cos(7/'2)---cos(7/'2'--i)v,sin(0i + Ci)...cos(0fe + Cfc) = sin(0i + 
Ci)...cos(0fc + Cfc) cos(-(/^fc+i)... cos(-?/'2'--i), where ipj = for each j > 1, 2 < 
< 2^ - 1, so that Tj cos(0i + (i) = for each j > 1 and / > 1, Tj sin(0i + 
Ci)...cos(0fc + Cfc) = for each j > A; and / > 1, where Tj^ = T]-^(Tj^) 
is the iterated composition for / > 1, / G N. Then T'^e""*^^'*'^-' gives with 
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such convention the same result as Sg e~'"^'*'''\ so one can use the symbohc 
notation Tje"'"*^'*'''^ = Q-uipM-ij'^i/'i) _ iq avoid misunderstanding we 

shall use Se^ and Tj in the sense of Formulas 12(3.1 — 3.7). 

It is worth to mention that instead of 12(3.7) also the formulas 

(1) exp(piii + ... +p„i„) = cos(0) +Msin(0) with := 0(p) := \p\-\- ...-\- 
p2]V2 and M = (piii + ... +pnin)l(t> for ^ 0, = 1; 

(2) 9exp(piii + .. + Pnin)/dpj = [-sin(0) + Mcos(0)]pj/0 + {(f)ij - 
Mpj)0~^sin(0) and d{pjtj + Cj)/^0 — 1 can be used. 

13. Theorem. Let f{t) be a function- original. Suppose that u{p, t; Q 
is given by 2(1,2,2.1) or i(8,8.1) over the Cayley-Dickson algebra Ar with 
2 < r < oo. Then a (super) derivative of an image is given by the following 
formula: 

(1) {^T-{f{t),u■,p■,0/^p).h^-T-{f{t)s^,u■,p■,C)ho 
-SeiJ""(/(t)si,it;p;C)/ii - - - Se^r'{f{t)sn,u;p;C)hn 
in the Ar spherical coordinates, or 

(1.1) {dT^{f{t),u;p;Oldp).h^-T^{f{t)s,,u;p;Oho 
-Se,J^^{f{t)t,,u;p;C)hi - ... - Se„T^{f{t)tn,u;p;C)hn 
in the Ar Cartesian coordinates for each h — /iq^o + ... + hnin £ Ar, where 
ho, hr, e R, 2''-^ <n<2'- -1, peWf. 

Proof. The inequalities ai(/) < Re{p) < a_i(/) are equivalent to the in- 
equalities ai{f{t)\t\) < Re{p) < a-i{f{t)\t\), since lim|(|_^+ooexp(— 6|i|)|i| = 
for each 6 > 0. An image T'^{f{t),u;p;C,) is a holomorphic function by p 
for ai(/) < Re{p) < a_i(/) by Theorem 4, also | /o°° e '*t'^dt\ < oo for each 
c > and n = 0,l,2,.... Thus it is possible to differentiate under the sign of 
the integral: 

(2) (a(/ /(t)exp(-«(p,t;C))dt)/ap)./i = 
^ {d{[ f{t)eM-<P,t-X))XuJt)/dp).h = 

v&{-l,l}^ -^^^ 

= / f{t){deM-<P,t;C))/dp).hdt. 

JR" 

Due to Formulas 12(3, 3.2) we get: 

(3) {dexp{-u{p, t; C))/dp).h = - ex.p{-u{p, t; C))si/io-Sei ex.p{-u{p, t; C))si^i- 
... - Se„ exp(-ii(p, t; C))s„hn 

in the Ar spherical coordinates, or 

(4) {dexp{-u{p, t; C))/dp).h = - ex.p{-u{p, t; C))si/io-Sei ex.p{-u{p, t; C))tihi- 
... - Se„ exp(-w(p, t; C))tnhn 

in the Ar Cartesian coordinates. 

Thus from Formulas (2, 3) we deduce Formula (1). 
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14. Theorem. If f{t) is a function-original, then 

(1) J^^ifit - r),u;p; C) = J^{f{t),u;p; C+ <p,r]) for either 
{{) u{p, t; C) = Posi + M{p, t; () + Co or 

(ii) u{p, t; C) =< p,t) + ( over Ar with 2 < r < oo in a domain p e Wf, 
where r e R", 2"^-^ < n < 2^ - 1, 

(2) < p,t] — poSi + piSiii + ... + PnSnin with Sj — Sj{n;T) for each j in 
the first (i) and < p, r] =< p, r) in the second (ii) case (see also Formulas 
1{8), 2(1,2,2.1)). 

Proof. For p in the domain Re{p) > Oi the identities are satisfied: 

(3) J'''{{fXm,..J{t-T),u;p;C)= H ... T f{t - T)e--^''-'iUt 

>'C/i,...,i 

due to Formulas 1(7, 8) and 2(1, 2, 2.1, 4), since poSi{n; t) + (o — poSi{n; ^) + 
Co+PoSi{n;r) and pjtj+Q ^ PjCj + {Cj+Pj^j) ^^dPjSj{n;t)+Cj ^ PjSj{n;^) + 
iCj + Pj^ji'iT'j '^)) for each j — 1, 2^ — 1, where i = ^ + r. Symmetrically 
we get (2) for instead oi\J\^...^\. Naturally, that the multiparameter non- 
commutative Laplace integral for an original / can be considered as the sum 
of 2" integrals by the sub-domains C/„: 

(4) / f{t)^-^v{-n{p,t;Q)dt= ^ / f{t)^-:^{-u{p,t;Q)xuX^)dt. 

The summation by all possible v e {—1, 1}" gives Formula (1). 

15. Note. In view of the definition of the non-commutative transform 
J*" and u{p,t;() and Theorem 14 the term (lii + ... + C2'--i«2'--i has the 
natural interpretation as the initial phase of a retardation. 

16. Theorem. If f{t) is a function- original with values in Ar for 2 < 
r <oo, 2"-^ <n <2'' -1, b eR, then 

(1) :F-{e''^'^+-+'-^f{t),u;p; Q = :F-{f{t),u;p- b; Q 

for each a_i -|- & > Re{p) > ai-\-b, where u is given by i(8, 8.1) or 2{1, 2). 

Proof. In accordance with Expressions 1(8,8.1) and 2(1,2,2.1) one has 
■"(P) ^; C) ~ ^(^1 + ■■■ + tn) — ^{p — b, t; C). If a_i + 6 > Re{p) > ai-\-b, then 
the integral 

(2) J-'^(e^(*^+-+*")/(i)xf/.(t),x.;p;C)= / /(t)e^(*^+-+*")exp(-x.(p,i;C))di 

= / f(t)eM-u(p-b,t;0)dt^r'(f(t)xuM,u;p-b;0 
converges. Applying Decomposition 14(4) we deduce Formula (1). 
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17. Theorem. Let a function f{t) be a real valued original, 

^ip'X) — -^"(/(O) '^iP! C); where the function u(p,t;Q is given by i(8, 8.1) 
or ^(1,2,2.1). Let also C) ''snc? q(p) be locally analytic functions such 
that 

(1) J^^{g{t, r);u;p; C) = G{p; C) exp(-ii(5(p), r; C)) 

for u =< p,t) +C oru^ po{ti + ... + t„) + M{p, t; C) + Co; then 

(2) ^"(/r. git, r)f{r)dT; u;p; C) = G{p; C)Fiqip); C) 

for each p e Wg and q{p) e Wf, where 2 < r < oo, 2^~^ < n < 2'' — 1. 

Proof. If p e Wg and e W/, then in view of the Fubini's theorem and 
the theorem conditions a change of an integration order gives the equahties: 

/ (/ g{t,T)f{T)dT)Q^v{-u{p,t-X))dt 
= / (/ g{t,T)eM-<P.t-X))dt)f{T)dT 

JR" JR" 

= / G'(p;C)exp(-^g(p),r;C))/(r)dT 

JR" 

= G(p;C) / /(T)exp(-i.(g(p),T;C))o?T = G(p;C)F(g(p);C), 

>/R" 

since r e R" and the center of the algebra Ar is R. 

18. Theorem. If a function f{t)xui,...,i is original together with its 
derivative d''f{t)xui,...,i{t)/dsi...dsn ord"'f{t)xui,...,i{'t)/dti...dtn, where F^{p;Q 
is an image function of f{t)xui,...,i over the Cayley- Dickson algebra Ar with 

2 < r e N, 2^-1 < n < 2"^ - 1, for u ^ PoSi + M{p, t; Q + Co given by 
£(1,2,2.1), then 

n-l 

(1) lim {bo +PlSeJp2Se....p„Se„F„"(p; () + E ("l)"^ 

m=0 

n—m 

<n; l<h<...<lm.<n; lai^jp Va,/3 

iC""^(p^'^C)} = (-l)"+V(0)e-"(°'°^«, 

or 

(1.1) lim {bo +PlSeJ[Po +P2Se.]...[po +PnSejF„"(p; C) + Y^i'^T 

m=0 

J2 \PO+Ph Se,i ] bo+Pj2 Se,-, ] ■ ■ ■ bo+Pj„-^ Se,_^ ] 

n—m <n; l<h<...<lm<n; lai^jp Va,/3 

^n-m(^(0.^)| = (_l)n+l/(0)e-"(O'O;^) 

for u{p,t; () given by i(8,8.1), where /(O) = lim^gc/j -^;t_^o/(^); P tends to 
the infinity inside the angle \Arg{p)\ < 7r/2 — 6 for some < 6 < it/2, 
l<j<2''- 1, p(') = J2]=o,j^(i)Pjij, (0 = ih, Im)- If the restriction 
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f{t)\tjj^=0,...,tj^=0;tk=oo\/k^{ji,...,jm} — l™te!7i,...,i;%i-^0,...,tj,„-s-0;tfc-^(X) Vfe^{ji,...,j^} f{t) 

exists for all 1 < ji < ... < < n, then 

n— 1 

(2) lim{[po+PlSeJp2Se....p„Se„F„-(p;C)+ E(-l)" 

m=0 

n — m 

<n; l<h<...<lm.<n; la^jp Vq:,/3 



n— 1 



— XK"!)"" X! /(^)|tji=0,...,tjv„=0;tfc=ooVfc^{il,...j„,}e 

m=0 l<j-i<...<jm<n 

in the Ar spherical coordinates or 



-«(o,o,c) 



n—l 



(2.1) lim{[po +PlSeJbo +P2Se,]...bo +PnSejF„"(P; + E i'^^ 

m=0 

E bo+Pji Se,, ] bo+Pj2 Se,2 ] - bo+Pj„-^ Se,_^ ] 

n — m 

<n; l<li<...<lm<n; laJ^j/S "^^,0 

n-l 

= E E /(Oltii=0,...,tj„=0;tfe=ooVfc^{il,...jVn}e""^°'°'^'' 

m=0 l<ji<...<jm<n 

m the Ar Cartesian coordinates, where p — >■ inside the same angle. 
Proof. In accordance with Theorem 12 the equahty follows: 

(3) J''^i{dfit)/ds,)xu,,...At),u;p-X) = bo5i,,+P,Se,]^"(/(t)xt/,,...,,(t),«(p,t;C),p;C) 

-r'-'''\mxu,_„u{p,t^;C);p;C) 

for u = u{p, t; C) = PqSi + M{p, t; () + Co iii the Ar spherical coordinates, or 

(3.1) r'((df(t)/dtj)xu^,..Jt),u;p;C) = \po+PjSe,]J^(mxu,,..Jt),u(p,t;0,P;0 

-J^-''''{f{t^)Xu,_„u{p,f;0;p-X) 

in the Ar Cartesian coordinates, since 

(3.2) dfit{s))/dsj = -df{t)/dtj^i+dfit)/dtj for each i > 2, dfit{s))/ds, = 
9/(t)/9ii, 

where p = po + ■■■ + P2'--i«2'--i e A, po, ■-,P2'--i e R, {io, ■-, «2'--i} 

are the generators of the Cayley-Dickson algebra Ar, Sn+i — for each Z > 1, 
the zero power 5°^ = / is the unit operator. For short we write / instead of 
/X!/i,...,i- Thus the limit exists: 

(4) j^^-'-''\f{f),u{p,P;cy,p;C)^ 

pco pco poo poo 

lini / dti... dtj-i dtj+i... (dtn)/(i) exp(-ii(p, t; C))- 

tj^+OJo Jo Jo Jo 
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Mention, that {...{{t^f)...y = (0, 0, tj, t„ : tj = 0) for every 1 < j < n, 
since = s^ — s^+i for each 1 < k < n. We apply these Formulas (3, 4) by in- 
duction j = l,...,n, T-^ <n< T-l,tod''f{t)/dsi...dsn,...,d''-^+^f{t)/dsj...dsn, 
...,df{t)/dsn instead of df{t)/dsj. 

From Note 8 [18j it follows, that in the Ar spherical coordinates 

p— >oo,|Arg(p)|<7r/2— 

also in the Ar Cartesian coordinates 

hm J^-{id-fit)/dh...dtn)xu,,..,„u;p;O = 0, 

p^QO,\Arg(p)\<TT/ 2—0 

which gives the first statement of this theorem, since u{p, 0, Q = u{0, t; Q = 
m(0, 0, C) and F°(j9(i--i); C) = /(0)e-"(°'°'^), while F„"(j9; () is defined for each 
Re{p) > 0. 

If the limit f{t'^-'^) exists, where t"^^-^ := (ti, ...,tn : tj = oo), then 

/■OO /"OO /"OO /"OO 

(5) lim / dti... / dtj-i / dtj+i... / {dtn)f{t)exp{-u{p,t;C)) 

=:r-~''<''>{f{t<^>),uip,t<^>;C);p;C). 

Certainly, {...{{t<^>)<^>)...)<^> = {ti,...,tn : ti = oo,...,tj = oo) for each 
1 < J < 'T'- Therefore, the limit exists: 



p^O,\Arg{p)\ 



lim ^ / (9"/(t)/5si...9s„)exp(-poSi-Co-M(p,t;C)) 
(a"/(t)/asi...as„)e-"(°'°'^)dt 



Ui 1 



Ujj=0,...,tj„=0;tfc=oo Vfc^{ji,...j„} 
m=0 l<ji<---<im<" 

= ni.^^i i[PO +PlSei]P2Se2.-PnSe„F„''(p;C) 

p— >0, 1 Arg(p) I <7r/2— (5 

71-1 



+ E(-ir 

m=0 

n — m 

<n; l<li<...<lm<n; laj^jp '^ct,l3 

i^r'"(/^;C) + (-i)"/(o)e-"(°'°'«}, 

from which the second statement of this theorem follows in the Ar spherical 
coordinates and analogously in the A,. Cartesian coordinates using Formula 
(3.1). 

19. Definitions. 

Let X and Y be two R linear normed spaces which are also left and right 
Ar modules, where 1 < r. Let Y be complete relative to its norm. We put 



23 



X®^ :— X ... X is the k times ordered tensor product over R of X. 
By Lq^k{X®^ ,Y) we denote a family of all continuous k times R poly-linear 
and Ar additive operators from X'^^ into Y. Then Lq^k{X'^^,Y) is also a 
normed R linear and left and right Ar module complete relative to its norm. 
In particular, Lq^i{X,Y) is denoted also by Lq{X,Y). 

We present X as the direct sum X — XqIq ® ... ® X2r-ii2r-i, where 
Xo,...,X2r-_i are pairwise isomorphic real normed spaces. If A e Lq{X,Y) 
and A{xh) — {Ax)h or A{hx) — h{Ax) for each x & Xq and h e Ar, then an 
operator A we call right or left Ar-lvcieox respectively. 

An R linear space of left (or right) k times Ar poly-linear operators is 
denoted by L^,fc(X®^y) (or Lr,k{X®^,Y) respectively). 

We consider a space of test function V :— I>(R",y) consisting of all 
infinite differentiable functions / : R" Y on R" with compact supports. 
A sequence of functions fn & V tends to zero, if all /„ are zero outside 
some compact subset K in the Euclidean space R", while on it for each 
A; = 0, 1, 2, ... the sequence {/^'^^ : n e N} converges to zero uniformly. Here 
as usually f'^'^^t) denotes the k-ih derivative of /, which is a times R poly- 
linear symmetric operator from (R")®*^ to Y, that is f^^\t).{hi, ...,hk) — 
/(*^)(i).(/io-(i), ^o-(fc)) £ ^ for each hi,...,hk G R" and every transposition 
a : {1, ...,k} — > {1, ■■■,k}, a is an element of the symmetric group S^, t e 
R". For convenience one puts /^^^ = /. In particular, f^^\t).{ej^, ■■■■>^jk) — 
d^f{t)/dtj^...dtj^ for all 1 < ji,...,jk < n, where ej = (0, 0, 1, 0, 0) e 
R" with 1 on the j-th place. 

Such convergence in T> defines closed subsets in this space T>, their com- 
plements by the definition are open, that gives the topology on T>. The space 
"D is R linear and right and left Ar module. 

By a generahzed function of class V :— [I>(R", Y)]' is called a continuous 
R-hnear ^^-additive function g : T> ^ Ar. The set of all such functionals 
is denoted by V. That is, g is continuous, if for each sequence /n £ "P, 
converging to zero, a sequence of numbers g{fn) —'■ [g, fn) £ Ar converges to 
zero for n tending to the infinity. 

A generahzed function g is zero on an open subset V in R", if [g', /) =0 for 
each f &T> equal to zero outside V. By a support of a generalized function 
g is called the family, denoted by supp{g), of all points t e R" such that 
in each neighborhood of each point t e supp{g) the functional g is different 
from zero. The addition of generalized functions g,his given by the formula: 

(1) [g + hj) := [gj) + [hj). 

The multiplication g eD' on an infinite differentiable function w is given 
by the equality: 
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(2) [gw, f) — [g, wf) either for w : R" — )■ Ar and each test function f eV 
with a real image /(R") C R, where R is embedded into Y; or w : R" — > R 
and / : R" ^ y. 

A generahzed function g' prescribed by the equation: 

(3) [g', f) :— —[g, /') is called a derivative g' of a generalized function g, 
where /' e P(R", L,(R", y)), g' e [P(R", L,(R", y))]'. 

Another space B :— B(R", Y) of test functions consists of all infinite dif- 
ferentiable functions / : R" — >■ Y such that the limit limui^^+oo 
exists for each m — 0,1, 2, j = 0, 1, 2, .... A sequence /„ e B is called con- 
verging to zero, if the sequence 1^1"*/^-'^^) converges to zero uniformly on 
R" \ S(R", 0, R) for each m,j = 0, 1, 2, ... and each < R < +oo, where 
B{Z,z,R) :— {y E Z : p{y,z) < R} denotes a ball with center at z of 
radius i? in a metric space Z with a metric p. The family of all R-linear and 

additive functionals on B is denoted by B'. 

In particular we can take X — A^, Y — A^ with 1 < a, /3 e Z. Anal- 
ogously spaces V(U,Y), [V(U,Y)]', B(U,Y) and [B(U,Y)]' are defined for 
domains U in R", for example, U — Uy (see also §1). 

A generalized function f e B' we call a generalized original, if there exist 
real numbers Oi < a_i such that for each oi < w^i, Wi, ty_n, < ci-i the 
generalized function 

(4) f(t)exp(-(qy,t))xu, is in [B(Uy,Y)y for all v = {vi,...,Vn), Vj E 
{— 1, 1} for every j — 1, ...,n for each t e R" with tjVj > for each j — 
1, ...,n, where = {viWy^i, 

By an image of such original we call a function 

(5) J""{f,u;p;() :— [/, exp(— C))) of the variable p e Ar with the 
parameter & Ar, defined in the domain Wf — {p & Ar '■ Oi < Re{p) < a_i} 
by the following rule. For a given p e Wf choose Oi < Wi,...,Wn < Re{p) < 
W-i, W-n < a_i, then 

(6) [/,exp(-'u(p,i;C)) := E^[/exp(-(g^, i)), exp{-['u(p, C)-(?^, i)]}X[/J, 
since exp{— t; C,) — {q^, t)]} e B{Uy, Y), where in each term 

[f exjp{-{qy,t)),exp{-[uip,t;C) - {qv,t)]}xuj the generalized function be- 
longs to [B{U-t,, Y)]' by Condition (4), while the sum in (6) is by all admissible 
vectors v e {—1, 1}". 

20. Note and Examples. Evidently the transform J""{f,u;p;Q does 
not depend on a choice of Wi, ty_„, ty„}, since 
[few(-((lv,t),ex.p(-[u(p,t;C) - {qv,t)])xu,) = 
[/exp(-(g^,i) - (by,t)),exp(-[u(p,t;C) - (qy,t) - (by,t)])xu,) 
for each b e R" such that oi < wj + bj < Re{p) < W-j -\- b^j < a_i for each 
j — l,...,n, because exp(— (6^, t)) e R. At the same time the real field R is 
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the center of the Cayley-Dickson algebra Ar, where 2 < r e N. 
Let 5 be the Dirac delta function, defined by the equation 



(DF) [S{t), (f){t)) := 0(0) for each (j) e B. Then 

(1) J^{5(^\t-r),u;p; C) = Eve{-i,i}4S^'Kt-r) exp(-(?,, f)), exp(-[ti(p, t; ()- 

{qv,t)])xuj 

^{-iy[dieM-Hp,t;cmt=r, 



since it is possible to take — oo < oi < < a_i < oo and Wk — for each k e 
{— 1, 1, —2, 2, — n, n}, where r e R" is the parameter, df :— d^^^/dfi ...dfi . 
In particular, for j' = we have 



(3) J^^{m{t)/^st..^si",u■,p■,0- 
^o<k,<n {j:yo'~''{PlSe,)''{P2Se,y'...{PnSeJ- exp(-Co - M(p, 0; ()) 



in the spherical coordinates, or 

(3.1) T-{d\^\5{t)/dt{\..dtt,u;p;0^ 

(PO +PlSeiyMPO +P2Se2y'...(po +PnSeJ^" exp(-u(p, 0; ()) 

in the Ar Cartesian coordinates, where ji + ■■■ + jn — ki, ji, jn are 
nonnegative integers, 2^~^ < n < 2'^ — 1, := l\/[m\(l — m)\] denotes the 
binomial coefficient, 0! = 1, 1! = 1, 2! = 2; Z! = 1 • 2 • ... • Z for each I > 3, 

The transform of any generalized function / is the holomorphic 

function by p G Wf and hy ( e Ar, since the right side of Equation 19(5) is 
holomorphic by p in Wf and by C in view of Theorem 4. Equation 19(5) im- 
plies, that Theorems 11-13 are accomphshed also for generahzed functions. 

For oi = a_i the region of convergence reduces to the vertical hyperplane 
in Ar over R. For a_i < ai there is no any common domain of convergence 
and f{t) can not be transformed. 

21. Theorem. If f{t) is an original function on R", F"{p-X) is its 
image, d\^^f{t)/ds^...ds^^ or d^^\ f{t) /di^ ...dP^ is an original, \j\ ^ 
3n, < Ji, e Z, 2^-1 < n < 2'- - 1; then 



(^^'jyo^"'HPlSej'HP2SeJ^'^..(PnSeJ^"^(/(i), C) 

for u{p, t; C) := Po^i + M{p, t; () + Co given by 2{1, 2, 2.1), or 

(1.1) :F^{d\^\f{t)/dt{\..dti",u-p-c)^ 

{po+PlSe,y'{po+P2Se,y'...{PO+PnSeJ-J^{f{t),U;p; Q 




(1) J'"{d\^^f{t)/ds^^...dsi-,u;p;0- E 
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for u{p,t; () given by i(8,8.1) over the Cayley- Dicks on algebra Ar with 2 < 
r < oo. Domains, where Formulas (1,1.1) are true may be different from a 
domain of the multiparameter noncommutative transform for f , but they are 
satisfied in the domain ai < Re{p) < a_i, where 

a_i = mm(a_i(/),a_i(aH/(t)/90r-..9C") : 1^1 < < < j,V/); 

ai = max(ai(/),ai(9H/(t)/50p...50™") : |m| < |A;|,0 < m; < if 
oi < a_i, where (pj = Sj or (pj = tj for each j correspondingly. 

Proof. To each domain 11^ the domain U-v symmetrically corresponds. 
The number of different vectors v G { — 1, 1}" is even 2". Therefore, for 
u = Pot + Co + M{p, t; due to Theorem 12 the equality 

(2) / (9/(t)/9s,)e-"(P'*'«)(is = / (9/(t)/9s,)e-"(P'*'«)(it = 

JR" JR" 



^{dtn[f{t)e--''''-'^^r^- ^{df){ f{t)[de-^^''-'^^/ds,]ds,) 

is satisfied in the Ar spherical coordinates, since the absolute value of the 
Jacobian dt/d{P , Sj) is unit. Since for ai < Re{p) < a_i the first additive is 
zero, while the second integral converts with the help of Formulas 12(2, 2.1), 
Formula (1) follows for A; = 1: 

(3) J^^{df{t)/dsj,u;p;0=po5i,,J^^{f{t),u;p;0+p,Se^J^''{f{t),u;p;0. 

To accomplish the derivation we use Theorem 14 so that 

lim[^"(/(t),n;p;C)-^"(/(t-re,),M;p;C)]/r 

r— >U 

= lim[^"(/(t), p; Q - J'''{f{t),u- p; C + r(po + PiH + ■■■+ P,^,))]/r 

= hm / /(t)[e-"(^''*'^) - e-"(P'*'C+^(P«+Pi^i+-+P^-^^))]r-Mt, 

where Cj = (0, 0, 1, 0, .., 0) G R*^ with 1 on the j-th place. If the orig- 
inal d^^'^f{t)/ds^...ds^^ exists, then 9l™l/(t)/9s™^--<9s™" is continuous for 
< \m\ < \j\ — 1 with < < ji for each I = 1, ...,n, where f^ := /. The 
interchanging of limT-_j.o and /j^„ may change a domain of convergence, but 
in the indicated in the theorem domain ai < Re{p) < a_i, when it is non 
void. Formula (3) is valid. Applying Formula (3) in the Ar spherical coor- 
dinates by induction to (9l'"l/(t)/9s™\..9s™") : 1^1 < |j|,0 < < j;V/) 
with the corresponding order subordinated to d^^^f{t)/ds\^...ds{^, or in the 
Ar Cartesian coordinates using Formula 12(1.1) for the partial derivatives 
(9l™l/()f:)/(9t™^---f^C") ■ 1""^! ^ lil)0 < "^z < jiV/) with the corresponding or- 
der subordinated to d^^^f{t)/dt-i...dP^ we deduce Expressions (1) and (1.1) 
with the help of Statement 6 from §XVII.2.3 [30] about the difi^erentiation of 
an improper integral by a parameter and §2. 
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22. Remarks. For the entire Euclidean space R" Theorem 21 for 
df{t)/dsj gives only one or two additives on the right side of 21(1) in accor- 
dance with 21(3). 

Evidently Theorems 4, 11 and Proposition 10 are accomplished for JF^'^^oW'-'^wi^f^ ^^-p- Q 
also. 

Theorem 12 is satisfied for and any j e {j(l), j(A;)}, so 

that si — si{k; t) — tj{i) + ... -\-tj{k) for each 1 <l < Pm — ^ and = for 
each 1 < m ^ 0(1); ■■■)i(^)} (the same convention is in 13, 14, 17, 21, see 
also below). For J^'^'^W'-'^Hk) in Theorem 13 in Formula 13(1) it is natural 
to put tjn — and — for each 1 < m ^ 0(1); so that only 

(k + l) additives with ho, hj{i),...,hj{k) on the right side generally may remain. 
Theorems 14 and 17 and 21 modify for putting in 14(1) and 

17(1, 2) and 21(1) tj — and Tj — respectively for each j ^ {i(l), ...,j{k)}. 

To take into account boundary conditions for domains different from U^, 
for example, for bounded domains V in R" we consider a bounded noncom- 
mutative multiparameter transform 

(1) T^{f{t)xv, u;p; C) =: T^{f{t),u-p- C). 
For it evidently Theorems 4, 6-8, 11, 13, 14, 16, 17, Proposition 10 and 
Corollary 4.1 are satisfied as well taking specific originals / with supports in 
V. 

At first take domains W which are quadrants, that is canonical closed 
subsets affine diffeomorphic with = n"=i[aj,^j], where — oo < ^j < bj < 
oo, [a^, bj] :— {x e H : < x < bj} denotes the segment in R. This 
means that there exists a vector w e R" and a linear invertible mapping C 
on R" so that C{W) - w ^ Q. We put P'^ := (ii, ...,tj, ...,tn : tj = a^), 
P''^ :— {ti, ...,tj, ...,tn : tj — bj). Consider t — {ti, ...,tn) £ Q"- 

23. Theorem. Let f{t) be a function- original with a support by t vari- 
ables in Q" and zero outside Q" such that df{t)/dtj also satisfies Conditions 
1{1 — 4). Suppose that u(p,t;Q is given by ^(1,2,2.1) or i(8,8.1) over Ar 
with 2<r <oo, 2'-^ < n < 2^ - 1. Then 

(1) J'''mit)/dt,)xQn{t),u;p;0^ 

J^-'-'''\m')XQr^(t^''),u;p; C) - J^-'-''''\m')xQr^(t^''),u;p; C) 

+bo + E PkSe,]:F''{f{t)xQ4t),u;p; Q 
k=i 

in the Ar spherical coordinates, or 

(1.1) J-((a/(i)/9i,)xQn(i),x.;p;C) = 

:F--'-^^^'\fif^^)XQr.if''),u;p;C)-r^-^^^^^^^^^ 
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+\po+PjSe,]J^{f{t)xQAt),u;p;C) 

in the Ar Cartesian coordinates in a domain W C Ar', if 3j ~ — oo or 
bj — +00, then the addendum with P'^ or P'^ correspondingly is zero. 

Proof. Here the domain is bounded and / is almost everywhere 
continuous and satisfies Conditions 1(1 — 4), hence f(t)exp{—u{p,t;()) e 
L^CFV^, Xn,Ar) for each p e Ar, since exp(— m(p, t; ^)) is continuous and 

suppifit)) c g". 

Analogously to §12 the integration by parts gives 
(2) f\df{t)ldt^)eM-<P,t-0)dtj = f{t))eM-<P,t\0)\X^% 



[/(t)(9exp(-«(p,t;C))/9t,)]rft„ 
where t — (ii, Then the Fubini's theorem implies: 

(3) / {df{t)/dt,)eM-u{p,t-Q)dt = 

■■■ / / [/ {df{t)/dt,)eM-<P.t;0)dt,]dP 
= [ f f{P'') exp{-u{p, P''; C))dP]-[ f f{P'') exp(-«(p, C))dP] 

JteQ^, tj=bj JteQ^, tj=aj 

•L^ rb\ rbn 

+ bo + Z^^^feSe,] / ... / f{t)exp{-u{p,t;C))dt 

k=l *^3n 

in the Ar spherical coordinates or 

(3.1) / {df{t)/dt^)eM-<P,t;C))dt 

Jt£Q", tj=bj Jt£Q", tj=aj 

fbl j-bn 

+bo+PjSeJ / ... / f(t)exp{-u(p,t;C))dt 

Jai Jan 

in the Ar Cartesian coordinates, where as usually P = {ti, 0, tj+i, ...,tn), 
dp — dti...dtj_idtj+i...dtn. This gives Formulas (1, 1.1), where 

(4) J^-'-'''-\fiP'')xQ<t^''),uip, P''; C);p; C) = 

fbl fbj~i fbj+i fbn 

... / / /(i^'^)exp(-i.(p,i^'^C))o?^''' 

J a\ Jaj — i Jaj^i Jan 

is the non-commutative transform by P''^, 2''^^ < < 2*" — 1, dP''^ is the 
Lebcsguc volume element on R"^^. 

24. Theorem. // a function f{t)xQ"{t) is original together with its 
derivative d^^ f{t)xQ^{t) / dsi...dsn or d^^ f{t)xQ^{t) / dti...dtn, where F^'{p-X) 
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is an image function of f{t)xQn{t) over the Cayley- Dickson algebra Ar with 
2 < r e N, 2^-1 <n<2'' -1, for the function u{p, t; () given by 2{l, 2, 2.1) 
or i(8, 8.1), Q" = n"=i[0, bj], bj > for each j, then 

(1) lim{[po+PlSeJp2Se,...p„Se„F„"(p;C)+ E(-l)"^ 

m=0 

E \PoSl,n + Ph ]Ph Se,-2 ■ ■ -Pjn-m Se^„_^ 

n—m 

<n; l<li<...<lm<n; laT^jp Vo,/3 

iC~"^(p^^^C)} = (-l)-+V(0)e-"(°'°^«) 
m ^/le spherical coordinates, or 

(1.1) lim {[po +PlSeJbo +P2Se.]...[po +PnSejF:(p; C) + E ("l)"^ 

m=0 

E [PO+Pji Se,.^ ] [P0+Pj2 Se,-, ] - [PO+Pj„_^ Se^_^ ] 

n—m 

m t/ie Cartesian coordinates, where /(O) = YmvteQ^, t^of{t), p tends to 
the infinity inside the angle \Arg{p)\ < tt/2 — 6 for some < 6 < 7r/2. 

Proof. In accordance with Theorem 23 we have Equahties 23(1,1.1). 
Therefore we infer that 

(2) T--'-''''\f{t^''')xQAt'^'),u{p, C);p; C) = 

hm / dti... I dtj-i / dtj+i... / (dtn)f(t)ex.p(-u{p,t;C)), 

where /3j-i = aj = 0, /3,-2 = 6^ > 0, A; = 1,2. Mention, that (...((t^''i)2''2)...)J^ = 
(t : ti = Pi^i^,...,tj = Pj^ij) for every 1 < j < n. Analogously to §12 we 
apply Formula (2) by induction j = 1, 2''^^ < 2'' — 1, to 

9-/(t(s))/asi...a9„,...,a--^+7(t(s))/9s,...as„,...,a/(t(s))/9s„ 

instead of df{t{s))/dsj, Sj = Sj{n;t) as in §2, or applying to the partial 
derivatives 

d^f{t)/dt^...dtn,...,d^-'^'f{t)/dt,...dtn,...,df{t)/dt^ 

instead of df{t)/dtj correspondingly. If Sj > for some j > 1, then Si > 
for Q" and limp^oo e""^^'*^"'^^ = for such t«, where t = (ti,...,t„), (/) = 
ih, Q, \l\ = /i + ... + In, t(') = {tf\ 4'^ = a, for = 1 and = 6,- 
for Ij = 2, 1 < j < 2^ — 1. Therefore, 

E (-l)l'l/(i«)e-«(P'*''';« = (-l)"/(0)e-"(°'°;^), 

^~^°°ije{l,2}; j=l,...,n 

since u{p, 0; C) = m(0, 0; C), where /((')) = limjeQ«;t^t(o /(t). 
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In accordance with Note 8 [IB] 



lim j^-{{d-f{t)/ds,...dsr.)xQ^{t),u{p,t- c);v,0 = Q 

in the Ar spherical coordinates and 

hm ^"((9'^/(t)/ati...9t„)xQ"W,w(p,t;C);p;C) = o 

p— >-oo,|yirg(p)|<7r/2— 

in the Ar Cartesian coordinates, which gives the statement of this theorem. 

25. Suppose that f{t)xQ"{t) is an original function, Q is its image, 

d^^^fit)XQ-it)/dt{\..dP- is an original, |j| = ji + ... + j„, < ji, j„ G Z, 
2*"^^ < < 2*^ — 1, — oo < afc < 6fc < oo for each A; = 1, n, (/) = (/i, 
/fc G {0, 1, 2}, = A for bounded Q". Let = G A : ai < i?e(p)} for 
bk = CO for some and finite for each k; W = {p ^ Ar Re{jp) < a_i} for 
3tk = — oo for some k and finite 6a: for each k\ W = {p E Ar '■ ai < Re{p) < 
a_i} when = — oo and bi = +oo for some k and /; t^'^ = (ti'', 

We put t^'^ = tk and = for = 0, tf'' = for Ik = 1, tf"* = bk for 
4 = 2, (g) = (gi, ...,g„), |g| = gi + ... + g„, 

ai = max(ai(/),ai(9l™l/(t)/9tr-^r) : 1^1 < |j|,0 < < j^Vfe), 

a_i = min(a_i(/),a_i(9l™l/(t)/9tr-5C") : H < \j\,0 < < 
jk V/c) if ai < a_i. 

If afc = — oo and 6^ = +oo for with a given k, then = 0. If either 
afc > — oo or bk < +oo for a marked k, then 4 G {0,1,2}. We also put 
hk = hk{l) = sign{lk) for each k, where sign{x) = —1 for x < 0, sign{0) = 0, 
sign{x) = 1 for x > 0, /i = h{l), \h\ = \hi\ + ... + \hn\, 

ilj) := {hsign{ji), ...,lnsign{jn)). 

Let the vector (/) enumerate faces dQ^i^ in dQ1_^ for \h{l)\ = A; > 1, so 
that dQl_^ = U|Moi=fe<5fo' ^ = for each (/) (m) (see also 

more detailed notations in §28). 

Let the shift operator be defined: 

T(m)F(p; := F{p; ( — {iirui + ... + z„m„)7r/2), also the operator 

(SO) S(„)F(p;C) :=S™^..S-"F(p;C), 
where (m) = (mi,...,m„) G [0, oo)" C R", S^^^^ = Sk{m) for each positive 
number < A; G R, So = / is the unit operator for (m) = (see also 
Formulas 12(3.1 — 3.7)). As usually let Ci = (1, 0, 0),...,e„ = (0, 0, 1) be 
the standard orthonormal basis in R" so that (m) = miCi + ... + m„e„. 

Theorem. Then 

(1) J'''{d\^\f{t)xQAt)/dt^...dti\u{p,t;Cy,p;0 = 
R^,;Rf^...R^,:^"(/(t)XQ.(t),«;p;C) 
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+ E 

l<\{lj)\; mk+qk+hk=jk;0<mk; 0<qk; hk=sign[h,jk); qk=0 foT lkjk=0, foT Cach k=l,...,n; (0e{0,l,2}'' 

for u{p, t; C) in the Ar spherical coordinates or the Ar Cartesian coordinates 
over the Cayley- Dickson algebra Ar with 2 < r < oo, where 

(1.1) Rei :=Po+PlSei; Re2 := Po + i^lSei + ]?2Se2,---, Re„ :=Po+PlSei + 

P2Se2 + ■■■ + Pn^e„ in the Ar sphcucal coordinates, while 

(1.2) := Po +PlSei; Re2 := PO +P2Se2;---; Re„ := Po + PnSe„ in the 

Ar Cartesian coordinates, i.e. Re^ = ^cjip) o^re operators depending on the 
parameter p. If tf'' = oo for some 1 < j < n, then the corresponding 
addendum on the right of (1) is zero. 

Proof. In view of Theorem 23 we get the equahty 

(2) / [(ai"^i+V(i)/9C'---5C\"'5C^'5Cr---^r)e""^'''*'^V^ = 
I (cii*^)[(aH/(i)/9i^i...ai;"")e-"^p'*;^)]|^t 

- / {dt^){t\d\"'\f{t)/dt'^\..dt';i^)[de-''^'^-^^^/dtk]dtk) 

is satisfied for < < for each k — l,...,n with |m| < On the 
other hand, for p e 1^ additives on the right of (2) convert with the help of 
Formula 23(1). Each term of the form 



/ (^ii('))[(al<'l/(i('))xaQ" (^^'))/aif...aee-"(^'*^^)l 



can be further transformed with the help of (2) by the considered variable 
tk only in the case Ik = 0. Applying Formula (2) by induction to partial 
derivatives d\^\ f / dt>^\ . .dt^ , d^^^-^' f /dti\..dt{- ,...,d^" f /dti",...,df /dtn as in 
§21 and using Theorem 14 and Remarks 22 we deduce (1). 

26. Theorem. Let f{t)xui i{t) be a function- original with values in 
Ar with 2<r <oo, T'-^ <n<2'' - I, u is given by 2{l, 2, 2.1) or 1{8, 8.1), 

(1) g{t) := / ... / f{x)dx, then 
Jo Jo 

(2) J-"(/xc7i At),'^;P;C) = Re^K,-Kr.J'''i9{t)xu, ,it),u;p;C) 

in the domain Re{p) > max(ai,0), where the operators Re^ are given by 
Formulas ^5(1.1,1.2). 

Proof. In view of Theorem 25 the equation 

(3) :F^{fxu, At):U;p;C)^ 

R,,R,,...Re^r\g{t),u;p;C) 
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+ E 

1<|«|; 0<mfc<l; mfe+?ifc=l; hk=sign(lk); for each k=l,...,n; qi=0,...,q„=0 

(-l)l«lR™^R-^..R-"^"-l'^«l((7(t«),M;p;C), 

is satisfied, since d"'g(t)/dti...dtn = {fXUi^...^i)(t), where ji = l,...,jn = 1, 
Ij = 1 for each j = l,...,n. Equation (3) is accomphshed in the same 
domain Re{p) > max(ai, 0), since g{0) = and g{t) also fulfills conditions of 
Definition 1, while ai{g) < max(ai(/),0) + b for each 6 > 0, where ai G R. 
On the other hand, g{t) is equal to zero on and outside in 

accordance with formula (1), hence all terms on the right side of Equation 
(3) with |/| > disappear and supp{g{t)) C f/i,...,i. Thus we get Equation 
(2). 

27. Theorem. Suppose that F^{p; () is animageJ^'''^^''"'^''{f{t)xui,...,iit),u]p 
of an original function f{t) for u given by 2{1, 2, 2.1) in the half space W := 
{p G Ar : Re{p) > fli} with 2 < r < oo, pi = 0,...,]9j„i = 0; Ci = 
n/2,...,(j^i = 7r/2 for each j > 2 in the Ar spherical coordinates or (i = 
0,...,(j-i = for each j >2 in the Ar Cartesian coordinates; 

(1) the integral f^l^^ F^{po + Z'.Qdz converges, where p = Po + Piii + 
... + pkik G Ar, Pj G R for each j = 0,...,2'' - 1, T'-^ < k < T - I, 
f/i,...,i := e R'' : ti > 0, > 0}. Let also 

(2) the function F^{p;Q be continuous by the variable p G Ar on the 
open domain W , moreover, for each w > ai there exist constants Cj > 
and > such that 

(3) \F''{p;C)\ < CJ exp{-e^\p\) for each p G Sn^ri), Sr := {z E Ar : Re{z) > 
w}, < R{n) < R{n + 1) for each n G N, lim„_^oo = cxo, where ai is 
fixed, ( = (oio + ... + Ckh £ Ar is marked, (j G R for each j = 0, ...,k. Then 

(4) n F\po + z- c)dz = 5_e,^'^*-'*H/(t)xc/.,...,,(t)/0, C), 

where pi = 0,...,pj^i = for each j > 2; d = 7r/2,...Xj-i = 7r/2 and 

~ '5j(A;;t) in the Ar spherical coordinates, while (i = 0,...,Ci-i = and 
C,j = tj in the Ar Cartesian coordinates correspondingly for each j > 1. 

Proof. Take a path of an integration belonging to the half space Re{p) > 
w for some constant w > ai. Then 

1 / f{t)exp{-u{p,t;C))dt\ <C [ exp{-{po~ai){ti + ...+tk))dt < oo 

converges, where C = const > 0, po > w. For tj > for each j = 1, k 
conditions of Lemma 2.23 [18] (that is of the noncommutative analog over 
Ar of Jordan's lemma) are satisfied. If tj — > oo, then Sj — )■ oo, since all 
ti,...,tk are non-negative. Up to a set of Lebesgue measure zero 
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we can consider that ti > 0,...,tk > 0. If Sj — > oo, then also Si — > oo. The 
converging integral can be written as the following limit: 

fooij r^'^j 

(5) / F%Q + z-X)dz= \mi F%o + z-X)QM-i^\A)dz 

for 1 < j < k, since the integral /-^5^[-^*^('U^ + -2; C)]'^-^ is absolutely converging 
and the limit limK^QCxp{—K,\z\) = 1 uniformly by z on each compact subset 
in Ar, where S is a purely imaginary marked Cayley-Dickson number with 
1 5* I = 1. Therefore, in the integral 

(6) / ' F\po + z-C)dz= \ f{t)[eM-<Po + z,t-X))]dt)dz 

the order of the integration can be changed in accordance with the Fubini's 
theorem applied componentwise to an integrand g — qqIq + ... + Qnin with 
fl'; e R for each I — 0, n: 

(7) / ' F\po + z;Odz^ dt{ \f{t)eM-y'iPo + z,t;0)dz) 

= / f{t){ \e-''^°+'''-^^\dz}dt. 

Generally, the condition pi = 0,...,pj_i = and Ci = 7i"/2,...,Cj-i = 't/2 
in the Ar spherical coordinates or (i = 0,...,(j-i = in the Ar Cartesian 
coordinates for each j' > 2 is essential for the convergence of such integral. 
We certainly have 

(8) cos{z;z^,+Cj)dz = [sin(^,e,+0)/0] 1^5^ = [- cos(^,e,+0+^/2)/0] (Q, 

Jpjtj 

and 

(9) sm{z*z^,+Cj)dzj = [- cos(^,^,+0)/e,] igj = [- sin(^,^,+0+7r/2)/e,] igj 

for each > and —oo<Pj<bj<oo and j = 1, ...,k. Applying Formulas 
(5 - 9) and 2(1, 2, 2.1) or 1(8, 8.1) and 12(3.1 - 3.7) we deduce that: 

n [F\po + z- C)]dz = S_e, / [f{t)/i^] exp{-ii(p, t- C)}dt 

= s_e,-F'='*^'-'*H/Wxc/i,...,iW/6-,^;p;C), 

where t = (ti, t^), Sj = tj + ... + tk for each 1 < j < k, = t^, = Sj in 
the Ar spherical coordinates or = tj in the Ar Cartesian coordinates. 

28. Application of the noncommutative multiparameter trans- 
form to partial differential equations. 

Consider a partial differential equation of the form: 
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(1) A[f]{t)=g{t), where 

(2) Am) :=E|,|<.a,(t)(al^l/(t)/9tf...at^"), 

aj(t) G Ak are continuous functions, where < k E Z, j = (ji,...,j„), 
|j| := ji + ... + jn, < jfe G Z, a is a natural order of a differential operator 
A, 2 < r, 2^'^ < n < 2^' — 1. Since Sk = Sk{n;t) = + ... + tn for each 
k = 1, ...,n, the operator A can be rewritten in s coordinates as 

(2.1) A[fMs)) := Eu\<Mt){d\^\f{t{s))/ds^^...dsi")- 
That is, there exists ^ for some j with |j| = a and hj = for \j\ > a, 
while a function X]j.|j|=a bj(t(s))si\..s{" is not zero identically on the corre- 
sponding domain V. We consider that 

(Dl) U is a canonical closed subset in the Euclidean space R"", that is 
U = cl{Int{U)), where Int{U) denotes the interior of U and cl{U) denotes 
the closure of U. 

Particularly, the entire space R" may also be taken. Under the linear 
mapping (ti, ...,tn) ^ (si, s„) the domain U transforms onto V. 

We consider a manifold W satisfying the following conditions {i — v). 

{{). The manifold W is continuous and piecewise C°, where C' denotes 
the family of / times continuously differentiable functions. This means by 
the definition that W as the manifold is of class C° fl Cj^^^. That is W is of 
class on open subsets Wqj in W and W \ (Uj Wqj) has a codimension 
not less than one in W. 

(ii). W = [jT=o where = [Jk Wo,k, Wj nWk = dl for each k ^ j, 
m = dim-RW, dimnWj = m — j, Wj+i C dWj. 

{Hi). Each Wj with j = 0, ...,m — 1 is an oriented C°-manifold, Wj is 
open in UfcLj W^fc- An orientation of W^j+i is consistent with that of dWj for 
each j = 0, 1, m — 2. For j > the set Wj is allowed to be void or non-void. 

(iv). A sequence W'^ of orientable manifolds embedded into R", 
a > 1, exists such that uniformly converges to W on each compact 
subset in R" relative to the metric dist. 

For two subsets B and E' in a metric space X with a metric p we put 

(3) dist{B, E) := maxjsup^g^ dist{{b}, E), supgg^; dist{B, {e})}, where 
dist{{b}, E) := infgg^; p{b, e), dist{B, {e}) := inf beB pip, e), b E B, e E E. 
Generally, dim-R.W = m < n. Let (ej^(x), eJi^(x)) be a basis in the 

tangent space T^W^ at a; G consistent with the orientation of W^., G N. 

We suppose that the sequence of orientation frames (ej(xfc), e^(xfc)) 
of at Xk converges to (ei(x), em(x)) for each x G Wq, where lim^ Xk = 
X G Wo, while ei{x),...,em{x) are linearly independent vectors in R*^. 

(f). Let a sequence of Riemann volume elements on W^ (see §XIII.2 
[3U] ) induce a limit volume element A on W, that is, X{BnW) = \imk^oo{B n 
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W'') for each compact canonical closed subset B in R", consequently, X(W\ 
Wq) = 0. We shall consider surface integrals of the second kind, i.e. by the 
oriented surface W (see (iv)), where each Wj, j = 0, ...,m — l is oriented (see 
also §XIII.2.5 [30J). 

Recall, that a subset V in R'^ is called convex, if from a,b & V it follows 
that ea + (1 - e)b G V for each e G [0, 1]. 

(vi) . Let a vector w G Int{U) exist so that U — w is convex in R" and 
let dU be connected. Suppose that a boundary dU of U satisfies Conditions 
{i — v) and 

(vii) let the orientations of dU^ and U'' be consistent for each G N (see 
Proposition 2 and Definition 3 |30]). 

Particularly, the Riemann volume element on dU'' is consistent with 
the Lebesgue measure on U'^ induced from R" for each k. This induces the 
measure A on dU as in (v). 

Also the boundary conditions are imposed: 

(4) f{t)\9u = fo{t'), {d\^\f{t)/dsf...ds1r)\9u = Mt') for |g| < a - 1, 
where s = (si, s„) G R", (g) = (gi, g„), |g| = gi + ... + g™, < G Z 
for each k, t E dU is denoted by t', /o, /(g) are given functions. Generally 
these conditions may be excessive, so one uses some of them or their linear 
combinations (see (5.1) below). Frequently, the boundary conditions 

(5) f{t)\au = fo{t'), {d'f{t)/diy%u = fi{t') for 1 </< a - 1 are also 
used, where v denotes a real variable along a unit external normal to the 
boundary dU at a point t' G dU^. Using partial differentiation in local 
coordinates on dU and (5) one can calculate in principle all other boundary 
conditions in (4) almost everywhere on dU . 

Suppose that a domain Ui and its boundary dUi satisfy Conditions (-D1, i— 
vii) and gi = gxui is an original on R" with its support in f/i. Then any 
original g on R" gives the original gxu2 ='■ 92 on R", where U2 = R" \ Ui. 
Therefore, gi + (72 is the original on R*^, when gi and g2 are two originals 
with their supports contained in Ui and U2 correspondingly. Take now new 
domain U satisfying Conditions (-D1, i — vii) and {D2 — D5): 

{D2) UdUi and dU C dUi; 

{D3) if a straight line ^ containing a point wi (see (vi)) intersects dUi 
at two points yi and y2, then only one point either yi or y2 belongs to dU, 
where Wi E Ui, U — Wi and Ui — Wi are convex; if ^ intersects dUi only at 
one point, then it intersects dU at the same point. That is, 

(DA) any straight line ^ through the point Wi either does not intersect 
dU or intersects the boundary dU only at one point. 

Take now g with supp{g) C U, then supp{gxui) C Ui. Therefore, any 
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problem (1) on Ui can be considered as the restriction of the problem (1) 
defined on [/, satisfying {Dl — DA, i — vii). Any solution / of (1) on U with 
the boundary conditions on dU gives the solution as the restriction f\ui on 
Ui with the boundary conditions on dU. 

Henceforward, we suppose that the domain U satisfies Conditions (Dl, DA, i— 
vii), which are rather mild and natural. In particular, for this means 
that either = — oo or bk — +oo for each k. Another example is: Ui is 
a ball in R" with the center at zero, U — Ui Li (R" \ C/i,...,i), wi — 0; or 
[/ = C/i U e R" : tn > -e} with a marked number < e < 1/2. But 
subsets dU(^i) in dU can also be specified, if the boundary conditions demand 
it. 

The complex field has the natural realization by 2 x 2 real matrices so 
that i = o) , = — ( i) ■ '^^^ quaternion skew field, as it is well-known, 
can be realized with the help of 2 x 2 complex matrices with the generators 

^ = (-1 o)' ^ = (o -i)' ^ = tro)' or equivalently by 4 x 4 real 

matrices. Considering matrices with entries in the Cayley-Dickson algebra Av 
one gets the complexified or quaternionified Cayley-Dickson algebras {Av)c 
or {Av)h with elements z — al + bi or z — al + bJ + cK + eL, where 
a, b,c,e e Avi such that each a & Ay commutes with the generators i, /, J, 
K and L. 

When r — 2, f and g have values in ^2 = H and 2 < n < 4 and 
coefficients of differential operators belong to A2, then the multiparameter 
noncommutative transform operates with the associative case so that 

for each a e H. The left hnearity property T'^{af) — aT'^{f) for any a e 
^j,K,L is also accomplished for either operators with coefficients in R or 
Ci = /R®iR or ^j,K,L = /R© JR©i^R®-LR and / with values in Av with 
1 < n < 2^^ — 1; or vice versa / with values in Ci or 'iilj,K,L and coefficients a 
in Av but with 1 < n < 4. Thus all such variants of operator coefficients 
and values of functions / can be treated by the noncommutative transform. 
Henceforward, we suppose that these variants take place. 

We suppose that g{t) is an original function, that is satisfying Conditions 
1(1 — 4). Consider at first the case of constant coefficients on a quadrant 
domain Q". Let be oriented so that = — 00 and b^ — -\-oo for each 
k <n — k; either = —00 or b^ — -\-oo for each k > n — k, where < k < n 
is a marked integer number. If conditions of Theorem 25 are satisfied, then 

(6) J'''{A[m),u;p;C)= ^ a,{[Re,(p)r[Re.(p)p...[Re„(p)?"^'^(/(i)XQn(i), C) 
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+ E 

l<I('i)l; rnk+qk+hk=jk\ 0<mfe; a<q^; hk=sign{lkjk)\ gk=0 for lkjk=0, for each fc=l,...,n; (Z)g{0,1,2}" 

^J^{g{t)xQn{t),u;p;C) 
for u{p,t; in the Ar spherical or Ar Cartesian coordinates, where the op- 
erators Rgjip) are given by Formulas 25(1.1) or 25(1.2). Here (/) enumer- 
ates faces (9(5"/) in for \h{l)\ = A; > 1, so that dQl_i = [j\h(i)\=kQ^i)^ 
nSQ"^) = for each (/) 7^ (m) in accordance with §25 and the notation 
of this section. 

Therefore, Equation (6) shows that the boundary conditions are neces- 
sary: 

(9l'?l/(t«)/atf ...ae)|,Q„^ for IjI < a, \m > 1, a, ^ 0, q, = for 
kjk = 0, mk + qk + hk = jk, hk = sign{lkjk), k = 1, n, t^^^ e dQ'^^y But 
dimndQ'^ = n - 1 for ^Q" 7^ 0, consequently, (<9l^l/(t('))/(9tf ...at;^")|aQ«^ 
can be calculated if know (5l^'/(t('))/(9t^|i)...(9t^^^))|aQ«^ for = where 
j3 = {Pi, (3m), m = \h{l)\, a number j{k) corresponds to /^(fe) > 0, since 
Qk — for Ik — and > only for l^jk > and k > n — k. That is, 
are coordinates in R" along unit vectors orthogonal to dQ'^^y 

Take a sequence of sub-domains C U^^^ C U for each e N so 
that each U'^ = \J^i^ Q^i is the finite union of quadrants Q'^i, m{k) e N. 
We choose them so that each two different quadrants may intersect only by 
their borders, each U'' satisfies the same conditions as U and 

(7) limfc^oo dtst{U, U'') = 0. 

Therefore, Equation (6) can be written for more general domain U also. 

For U instead of wc get a face c'[/(/) instead of dQ'^i^ and local coordi- 
nates r^,(i),...,r^(m) orthogonal to dU(^i) instead of t^(i),...,t^(j„) (see Conditions 
(i — iii) above). 

Thus the sufficient boundary conditions arc: 

(5.1) {d\^\fm/dr%..3T'^-^^^^^^^ 

for = where m = \h{lj)\, \j\ < a, > 1, ^ 0, Qk = for 

hjk = 0,mk + qk + hk= jk, hk = sign{lkjk), < < jfe - 1 for A; > n - 
<Pi3,ii)i'^^''^) known functions on dU^j, t^^^ G dU{i). In the half-space tn > 
only 

(5.2) d^f{t)/dtl^=o 

are necessary for /3 = |g| < a and q as above. 

Depending on coefficients of the operator A and the domain U some 
boundary conditions may be dropped, when the corresponding terms vanish 
in Formula (6). For example, ii A = d'^/dtidt2, U = Ui^i, n = 2, then 
df /dulouo is not necessary, only the boundary condition f\Qu is sufficient. 
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U U — R", then no any boundary condition appears. Mention that 

(5.3)^0(/(a);i.;p;C) = /(a)e-"(^''«'O, 
which happens in (6), when a — t^'^ and \h{l)\ — n. 

Conditions in (5.1) are given on disjoint for different (l) submanifolds 
dU^) in dU and partial derivatives are along orthogonal to them coordinates 
in R", so they are correctly posed. 

In Ar spherical coordinates due to Corollary 4.1 Equation (6) with dif- 
ferent values of the parameter ( gives a system of linear equations relative to 
unknown functions S(^)J-""'(/(i), C), from which J^^{f{t),u;p-X) can be 
expressed through a family 

and polynomials of p, where Z denotes the ring of integer numbers, where the 
corresponding term jg zero when t^-^ — ±oo for some j. In the Ar 

Cartesian coordinates there are not so well periodicity properties generally, so 
the family may be infinite. This means that T"{f{t),u;p; () can be expressed 
in the form: 

(8) J'^{f{t),u;p;C) = J2Pim)iP)Sim)J'''{9{t),u;p;C) 

(m) 

where P(m,)(p) and Pj,{q),{i),(m){p) are quotients of polynomials of real variables 
Po,Pi, ■■■,Pn- The sum in (8) is finite in the Ar spherical coordinates and may 
be infinite in the Ar Cartesian coordinates. To the obtained Equation (8) wc 
apply the theorem about the inversion of the noncommutative multiparam- 
eter transform. Thus this gives an expression of / through (7 as a particular 
solution of the problem given by (1,2,5.1) and it is prescribed by Formulas 
6.1(1) and 8.1(1). 

For -F"(/; u; p; Q Conditions 8(1, 2) arc satisfied, since P(m)(p) and Pj.{q).{i).{m){p) 
are quotients of polynomials with real, complex or quaternion coefficients and 
real variables, also G" and terms on the right of (6) satisfy them. 

Thus we have demonstrated the theorem. 

28.1. Theorem. Suppose that J-'^{f;u;p;() given by the right side of 
(8) satisfies Conditions 5(3). Then Problem (1,2,5.1) has a solution in the 
class of original functions, when g and 4'fj^{i) are originals, or in the class of 
generalized functions, when g and (f)j3^(i) are generalized functions. 

Mention, that a general solution of (1,2) is the sum of its particular 
solution and a general solution of the homogeneous problem Af = 0. If 
0/3,(0 = + 9 = 9i + 92, f = /i + /2, Afj = gj and on dUj 



39 



satisfies (5.1) with j = 1, 2, then Af = g and / on dU satisfies 

Conditions (5.1) with 0/3,(;). 

28.2. Example. We take the partial differential operator of the second 
order 

n n 

A= ^h,md'^ / dThdTm + <^hd/dTh + 

h,m=l h=l 

where the quadratic form a(r) := J2h.m^h,mThTm is non-degenerate and is 
not always negative, because otherwise we can consider —A. Suppose that 
SLh,m = SLm,h £ R-, Cih,Th € R for each h,m = l,...,n, u G ^3. Then we 
reduce this form a(r) by an invertible R linear operator C to the sum of 
squares. Thus 

n n 

(9) A=J2^hdydtl + J2M/dh + u;, 

h=l h=l 

where (ti, t„) = (ri, r„)C with real and f3h for each h. If coefficients 
of A are constant, using a multiplier of the type exp{J2h ^hSh) it is possible 
to reduce this equation to the case so that if a.h 7^ 0, then (3h = (see §3, 
Chapter 4 in [26]). Then we can simplify the operator with the help of a linear 
transformation of coordinates and consider that only /3„ may be non-zero if 
a„ = 0. For A with constant coefficients as it is well-known from algebra one 
can choose a constant invertible real matrix {ch,m)h,m=i,...,k corresponding to 
C so that a/i = 1 for < k+ and a/^ = — 1 for /i > A;_|_, where < A;+ < n. 
For k+ = n and /3 = the operator is elliptic, for A;+ = ri — 1 with a„ = 
and /3n ^ the operator is parabolic, for < A;_|_ < n and (3 = the operator 
is hyperbolic. Then Equation (6) simplifies: 

n 

(10) J-"(A[/](t),n;p;C) = E a4[Re.(p)]'^'^(/(t)XQ"(t),«;p; 

h=l 

+ E (-l)^M-^""'W(i^'^)xao'^^,(t«)/at.,u;p;C) 

«he{l,2};(0=ihe,, 

+ [Re,(p)]^"-^(/(t«)XaQ'^',(t«),n;p;C)]} 
+/3„{^"-^^*"''(/(r'2)xaQ^^,„ (r-^), u; p; () -^"-i^*""(/(r'i)xaQj„ it-''), u- p; () 
+ [Re„(p)]J-"(/(t)xQn(t),n;p;C)}+u;^"(/(t)XQ"(t),u;p;C) = -^"(^7(t),«;p;C) 

in the Ar spherical or Ar Cartesian coordinates, where e/i = (0, 0, 1, 0, 0) G 
R" with 1 on the h-th place, Sq = / is the unit operator, the operators Re^ (p) 
are given by Formulas 25(1.1) or 25(1.2) respectively. 

We denote by Ss{x) the delta function of a continuous piecewise differen- 
tiable manifold 5* in R" satisfying conditions {i — vi) so that 

(A) / ri{x)5s{x)dx = / i]{y)X^{dy) 
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for a continuous integrable function rjlx) on R", where dim{S) — m < n, 
^m{dy) denotes a volume element on the m dimensional surface ^S" (see Con- 
dition {v) above) . Thus we can consider a non-commutative multiparameter 
transform on dU for an original f on U given by the formula: 

(11) J-a"f7'^*'(/(0xaf/(0,^;p;C) =-^"'*(/(^)<^af/(i),^^;p;C)- 

Therefore, terms hke J"""^ in (10) correspond to the boundary dQ'". They 
can be simplified: 

(12) /3n{-F"-^^*"''(/(r'2)xaQ5eJ^)'«;p;C)--^-'^*"'(/(^"'')xaQ?Jt),^ 

where (3{t') is a piccewise constant function on dQ""" equal to /3„ on the cor- 
responding faces of Q" orthogonal to given by condition: either tn = a„ 
or tn = bn] /3{t') = is zero otherwise. 

If = — oo or bk = +oo, then the corresponding term disappears. If R" 
embed into Ar with 2"^^^ < n < 2'' — 1 as Hii © ... © Ri„, then this induces 
the corresponding embedding of Q"^ or U into Ar- This permits to make 
further simplification: 

n 

(12.1) ^a,{ Y: (-i)'M[Re.(p)]-F^^-^(/(t«)xaQ5^(t('^«;p;C) 

h=l lh&{l,2}-{l)=lheh 

+^"-^(a/(t«)X9Q^,^(t«)M,,«;p;C)]} 
= J^Q^(a(t')(df(t')xdQ^^(t')/du),u(p, t'- C);p; C) 

+-^aQ»(P(^')/(OxaQ;^(0,«;p;C), 

where v = v{t') denotes a real coordinate along an external unit normal M{t') 
to Q{dU) at Q{t'), so that M{t') is a purely imaginary Cayley-Dickson num- 
ber, a(t') is a piecewise constant function equal to SLh for the corresponding 
t' in the face dQl,^ with > 0; P{t',p) := P{t') := K^p) for t' G dQ^,^, 
h = l,...,n, since sin(-?/; + vr) = —sm[ijj) and cos(-?/' + tt) = — cos('0) for 
each -0 e R. Certainly the operator-valued function P{t') has a piecewise 
continuous extension P{t) on Q"'. That is 

(13) J'Su\mfit')xau{t'),u{p,t';Cy,p;C) 

:= / e(^)/(^)<^ac/(i)exp{-ti(p,f;C)}dt 

for an integrable operator-valued function ^{t) so that [^(t)f(t)] is an original 
on U whenever this integral exists. For example, when ,^ is a linear combina- 
tion of shift operators S(^„i) with coefficients e(„j)(t,p) such that each ^(rn){t,p) 
as a function hj t E U for each p E W and f{t) are originals or / and g are 
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generalized functions. For two quadrants Qm,i and Qm,k intersecting by a 
common face T external normals to it for these quadrants have opposite di- 
rections. Thus the corresponding integrals in J^qqI^ ^ and J^q^^ restricted 
on T summands cancel in J-'a^ , .n v 

Using Conditions {iv — vii) and the sequence and quadrants out- 
lined above we get for a boundary problem on U instead of the following 
equation: 

n 

(14) r\A[f]{t),u-p-C) = {Ea.[Re,(p)]V"(/(t)xt/(t),«;p;C)}+ 

h=l 

where P{t',p) := P(t') := El=i Sih[^h{p)]{du / dtt) for each t' G dUo (see also 
Stokes' formula in §XIII.3.4 [50] and Formulas (14.2,14.3) below). Partic- 
ularly, for the quadrant domain we have a{t) = a.^ for t G dQ^^^^ with 
Ih > 0, (3{t) = Pn for t e c?Qr„e„ with /„ > and zero otherwise. 
The boundary conditions are: 

(14.1) f{t)\au, = mauo, {df{t)/du)\9Uo = Mt)\ouo- 

The functions a(t) and (3{t) can be calculated from {a^ : h} and /3„ almost 
everywhere on dU with the help of change of variables from (ti,...,t„) to 
(?/i, ?/„_i, Un), where (yi, ?/„) are local coordinates in dUo in a neighbor- 
hood of a point t' G Sf/o; 1/n = ^, since 9[/o is of class C^. Consider the 
differential form J2h=ii^^T~''^hdti A ... A dt^ A ... A dtn = adyi A ... A dyn-i 
and its external product with du = Ylh^iidu / dth)dth, then 

(14.2) a{t)\9u, = ELi Mdi^/dth)\dUo and 

(14.3) (3{t)\dUo = /3nXu.„uU2e„idu/dt^)\9Uo- 
This is sufficient for the calculation of J^gjj^. 

28.3. Inversion procedure in the Ar spherical coordinates. 

When boundary conditions 28(5.1) are specified, this equation 28(6) can 
be resolved relative to J^"'{f(t)xu(t),u{p,t;();p;(), particularly, for Equa- 
tions 28.2(14, 14.1) also. The operators Sg^. and Tj of §12 have the period- 
icity properties: SifF{p;0 = S':F{p;C) and T}+''F{p;0 = T^F{p;0 , 
Sl+^F{p- C) = -S^e,F{p; C) and T^^''F{p; () = -T^F{j)- Q for each positive 
integer number k and 1 < j < 2*' — 1. We put 

(6.1) Y,{p- C) := (S^^ - S%jF{p- C) for any 1 < j < 2^ - 2, 

(6.2) F2r_i(p;C) := Sg2^_^F(j9; C). Then from (6) we get the following 
equations: 

(6.3) a,{h+PiTipbo+Piri+p2T2p 
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...\po+PlTi + ... + PnTnY"}\p,=0 Vb>w F^iV, C) = {" E 

lil<a 

E 

1<I(0')I; rnu+qk+hk=jk; 0<mfc; 0<qk; hk=sign{lkjk); gk=0 for lkjk=0, for each fc=l,...,n; (i)e{0,l,2}' 

for each w = l,...,n, where F{p;() = T"'{f{t)xQ^^{t),u;p;() and G{p;() = 
J-'"'{g{t)xQ"{t),u;p; (). These equations are resolved for each w = 1, ...,n as 
it is indicated below. Taking the sum one gets the result 

(6.4) F(p;C) = Fi(p;C) + ... + F4p;C), 
since {[T^lfiSt^ - St^J] + S^^,_ Je-"(^'*^0 = Si^e-^^^'^-'O = e-«(^''*^C). 

The analogous procedure is for Equation (14) with the domain U instead 
of Q^ 

From Equation (6.3) or (14) we get the linear equation: 

(15) E V'(z)a;(z) = 0, 
(0 

where is the known function and depends on the parameter (, are 
known coefficients depending on p, X{i) are indeterminates and may depend 
on li = 0, 1 for h — 1, so that X(i)^2ei = = 0, 1,2,3 for h > 1, 

where a;(;)+4e^ = xq^ for each /i > 1 in accordance with Corollary 4.1, (/) — 

(/l, In)- 

Acting on both sides of (6.3) or (14) with the shift operators T(m) (see 
Formula 25(5*0)), where mi = 0,1, rrih = 0,1,2,3 for each /i > 1, we get 
from (15) a system of 2^+^'^'^""^) linear equations with the known functions 

4>{m) ■= T(^m)4> instead of 0, (/!)(o) = 4>: 

(15.1) E{i)i'ii)T(m)X{i) = (f){m) for each (m). 

Each such shift of ( left coefficients ■0(0 intact and X(^i)+(^m) = {—^Y'^H') 
with I' I = /i + mi {mod 2), 1'^ = Ih + iTih {mod 4) for each h > 1, where 
77 = 1 for /i + mi — /'i = 2, ?7 = 2 otherwise. This system can be reduced, 
when a minimal additive group Q := {(/) : h {mod 2), Ij {mod 4) V2 < j < 
k; generated by all {I) with non-zero coefficients in Equation (15)} is a 
proper subgroup of g2 x g^^^, where gh := Z/ {hZ) denotes the finite additive 
group for < /i e Z. Generally the obtained system is non-degenerate for 
A„+i almost all p = (po, •••,Pn) £ R""*"^ or in W, where A„+i denotes the 
Lebesgue measure on the real space R"+^. 

We consider the non-degenerate operator A with real, complex Ci or 
quaternion Hjx,l coefficients. Certainly in the real and complex cases at 
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each point p, where its determinate A = A(p) is non-zero, a solution can be 
found by the Cramer's rule. 

Generally, the system can be solved by the following algorithm. We can 
group variables by h, l2,---,lk- For a given I2, Ih and Zi = 0, 1 subtracting 
all other terms from both sides of (15) after an action of TJ^^) with mi = 0, 1 
and marked nih for each h > 1 we get the system of the form 

(16) axi + /3x2 = bi, 
—/3xi + ax2 = b2, 

which generally has a unique solution for A„+i almost all p: 

(17) X, = {a{a^ + P')-^)h, - {(3{a' + (3')-^)h2y, X2 = {a{a^ + P'r')h2 + 

where bi, 62 G Ar, for a given set (m2, ...m^). 

When Ifi are specified for each 1 < h < k with h ho, where 1 < ho < k, 
then the system is of the type: 

(18) axi + bx2 + CX3 + dx4 — bi, 
dxi + ax2 + bxs + 0x4 — b2, 

cxi + dx2 + axs + hxi — ba, 

hxi + CX2 + dx^ + = b4, 
where a,b,c,d e R or Ci or 11j^k,l, while bi, b2, bs, b4 e Ar- In the latter 
case of tlj^K,L it can be solved by the Gauss' exclusion algorithm. In the first 
two cases of R or Ci the solution is: 

(19) Xj = Aj/A, where 

A = a^i - dC2 + 0^3 - Ki, 

Ai = bi^i - \i2^2 + bs^s - b4^4, 
A2 = -bi^4 + ^2^1 — b3^2 + b4^3, 
A3 = bi^3 - \i2iA + b3^i - b4^2, 
A4 = -bi^2 + ^2^3 - ^Z^A + b4^i, 

^1 = + + cd? - ac^ - 2abd, 
^2 = a% + h(? + d^ - h'^d - 2acd, 
^3 = + c^ + ad^ - a^c - 2bcd, 
^4 = a'^d + b^ + c^d - hdP - 2abc. 

Thus on each step either two or four indeterminates are calculated and 
substituted into the initial linear algebraic system that gives new linear alge- 
braic system with a number of indeterminates less on two or four respectively. 
May be pairwise resolution on each step is simpler, because the denominator 
of the type (a^ -|- ^^) should be A2'- almost everywhere by p e positive (see 
(6), (14) above). This algorithm acts analogously to the Gauss' algorithm. 
Finally the last two or four indeterminates remain and they are found with 
the help of Formulas either (17) or (19) respectively. When for a marked h 
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in (6) or (14) all Ih — {mod 2) (remains only xi for h — 1, or remain xi 
and ^3 for /i > 1) or for some h > 1 all Ih — {mod 4) (remains only xi) a 
system of linear equations as in (15, 15.1) simplifies. 

Thus a solution of the type prescribed by (8) generally Xn+i almost ev- 
erywhere by p e 1^ exists, where 1^ is a domain W — {p e Ar ■ ai < 
Re{p) < a_i, pj — yj > n} of convergence of the noncommutative multi- 
parameter transform, when it is non-void, 2^~^ < n < 2^ — 1, Re{p) — po, 

P^Poio + ■■■ +Pnin- 

This domain W is caused by properties of g and initial conditions on dU 
and by the domain U also. Generally U is worthwhile to choose with its 
interior Int{U) non-intersecting with a characteristic surface 0(xi, x„) = 
0, i.e. at each point x of it the condition is satisfied 

{cs) j:\j\=aMt{x)){d<f>/dx,yK..{dci>/dx„y-^o 

and at least one of the partial derivatives {d(f)/dxk) 7^ is non-zero. 

In particular, the boundary problem may be with the right side g — qf in 
(2, 2.1, 14), where ^ is a real or complex Ci or quaternion 11j,k,l multiplier, 
when boundary conditions are non-trivial. In the space either Z>(R", Ar) 
or S(R",^r) (see §19) a partial differential problem simphfies, because all 
boundary terms disappear. If / G B(R",^r), then {p E Ar ■ Re{p) > 0} C 
Wf. For / e ©(R", A) certainly Wf = A (see also §9). 

28.4. Examples. Take partial differential equations of the fourth or- 
der. In this subsection the noncommutative multiparameter transforms in 
Ar spherical coordinates are considered. For 

(20) A = d^/dsf + E"=2 Ijd^/ds^j 
with constants 7^ e Hj,k,l\{0} on the space either D(R", Ar) or B(R", Ar), 
where n > 2, Equation (6) takes the form: 

(21) J^-{A[f]{t),u;p;C)^ 

n 

{P0(p^+3(piSej')+E7i(P,Se,)n-^"(/W,«;P;C)+Pl(3p^+(piSej')Se,.F"(/(t),«;p;C) 

^J^{g{t),u-p-C) 

due to Corollary 4.1. In accordance with (16, 17) we get: 

(22) C) = {a{a^ + P^)-')G^{p; Q - (/?(«' + P^)-')T,G^{p; Q) for 

each w = 1, n, 

where = a = [pq{pI - Spj) + Ej=2ljPj]\pi=o vb>w, Pw = P = Pii'ipl - 
Pi)lp6=o V6>u;- From Theorem 6, Corollary 6.1 and Remarks 24 we infer that: 

(23) f{t) = {2n)-^[ F{a+p;C)exp{u{p,t;C)}dpi...dpn 
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supposing that the conditions of Theorem 6 and Corollary 6.1 are satisfied, 
where F{p;C) ^ J^im,u;p;C). 

If on the space either T>(R'^, Ar) or B{B}, Ar) an operator is as follows: 

(24) A = OyOsldsl + E"=3 7^-9794, where 7,- e IIj,k,l \ {0}, where 
n > 3, then (6) reads as: 

(25) J-(A/(i),«;p;C) =pi(p^ + (piSej2)S^^J-"(/(t),«;p;C) 
+2poPiplS,,SlT^(f(t)),u;p; () + T.U^jip^S,,fT-{f{t)),u-p- C) 

^F^{g{t),u-p-C). 

If on the same spaces an operator is: 

(26) A = d^/dsidsl + E]=3ljd^/ds^j, where n > 3, then (6) takes the 
form: 

(27) T^{Af{t),u;p; C) = PopIS^, J-"(/(t), C)+PiPlSe,SlT^{m,u;p- ()+ 

j:Ulj{PjSe,)'J^{m,u;p;C)^J'''{g{t),u-p;C). 

To find T'^{f{t),u;p; Q in (25) or (27) after an action of suitable shift op- 
erators T'(o,2,o,...,o), ^(i,o,...,o) and T'(i,2,o,...,o) we get the system of hnear algebraic 
equations: 

(28) axi + bx3 + CX4 — bi, 
bxi + CX2 + ax3 = b2, 

0x2 — cxs + bx4 — ba, 

—cxi + bx2 + ax4 — b4 
with coefficients a, b and c, and Cayley-Dickson numbers on the right side 
bi, b4 e Ar, where Xi = F^(p; C),X2^ TiF^{p; C),Xs^ 7f F^(p; C), 2:4 = 
T,TiF^(p;0, bi = G^(p;C) = (J'''(9(t),u;p;0U ^2 = 7fG^(p;C), = 
riG^(,(p; C), b4 = TiTf G^(p; ()■ Coefficients are: a^^, = a = [E"=3 7jP^] Ip6=o V6>i(; e 
Hj,ft:,L, = 6 = pUPq-Pi) G R, = c = 2poPipi|pi,=o V6>«; e R for A given 
by (24); = a = [E"=3 TiP^] Ip6=o \/b>w e Hj^^.l, 6^ = 6 = PoPi|pi,=o \/b>w e 
R, = c = PiP2lpi,=o V6>w e R for A given by (26), w — 1, n. If a = the 
system reduces to two systems with two indeterminates {xi,X2) and (a;3,a;4) 
of the type described by (16) with solutions given by Formulas (17). It is 
seen that these coefficients are non-zero A„+i almost everywhere on R""*"^. 
Solving this system for a 7^ we get: 

(29) F^(p; C) = a-^bi - [a^ - 6^ + + 4i,2^2]-i^-i[(^2 _ ^2 ^ ^2)((^2 _ 

6^)bi a6b2 - 26cb3 + adOi) - 26c(26cbi - acb2 + (c^ - 6^)b3 + abh^)]. 
Finally Formula (23) provides the expression for / on the corresponding 
domain W for suitable known function g for which integrals converge. If 
7j > for each j, then a > for each p\-\- ... -\- p^^ > {). 

For (21, 24) on a bounded domain with given boundary conditions equa- 
tions will be of an analogous type with a term on the right T"'{g{t),u;p; Q 
minus boundary terms appearing in (6) in these particular cases. 
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For a partial differential equation 



(30) a(t„+i)yl/(ti, tn+i) + df{ti, ...,tn+i)/dtn+i = gih, tn+l) 

with octonion valued functions /, g, where ^4 is a partial differential operator 
by variables ti, i„ of the type given by (2, 2.1) with coefficients independent 
of ii, tn, it may be simpler the following procedure. If a domain V is not 
the entire Euclidean space R""*"^ we impose boundary conditions as above 
in (5.1). Make the noncommutative transform J^";*i> - >*" of both sides of 
Equation (30), so it takes the form: 

(31) a(i„+i) ^•-■*i'-'*"(A/(ii, tn+,),u;p; C)+dJ^'''-'-{f{h, tn+,),u;p; C)/9i„+i 

^J^''''-'''^{g{h,...,tn+i),u;p;C). 
In the particular case, when 

for each tn+i, p, t and (, with the help of (6, 8) one can deduce an expression 
of F"(p;C;W) - T^''''-''-(f(h,...,tn+i),u;p;0 through G^(p;C;tn+i) := 
J'"''*i' - '*"(g'(ti, tn+i),u;p; C) and boundary terms in the following form: 

(32) h{po,...,Pn]tn+l)F''{p]C]tn+l)+dF''{p]C]tn+l)/dtn+l = Q{po, ■■ Pn]tn+l) , 

where b(po, ■■■,Pn',tn+i) is a real mapping and Qipo, ■■■,Pn',tn+i) is an octo- 
nion valued function. The latter differential equation by tn+i has a solution 
analogously to the real case, since tn+i is the real variable, while R is the 
center of the Cayley-Dickson algebra Ar- Thus we infer: 

(33) F"(p; C; tn+i) = exp{- b(po, -.Pn; 

J To 

rt„+i rT 

{Co+[/ (5(po, -,Pn;T)exp{ / h{po,...,Pn]C)di}dT\), 

Jto Jto 
since the octonion algebra is alternative and each equation hx = c with non- 
zero b has the unique solution x = b^^c, where Co is an octonion constant 
which can be specified by an initial condition. More general partial differen- 
tial equations as (30), but with d'' f /dt''^^-^, / > 2, instead of df/dtn+i can be 
considered. Making the inverse transform (^^»*;*iv,<n^-i ^j^^ right side of 

(33) one gets the particular solution /. 

28.5. Integral kernel. We rewrite Equation 28(6) in the form: 

(34) A5J^(/xQn,ii;p;C) =^(^XQ",^i;p;C)- 
Ea. E 

lil<o; l<\{lj)\, 0<mk, 0<qk, hk=sign{lkjk), mk+gk+hk=jk; 9fc=0 for lkjk=0; VA;=l,...,n; (0e{0,l,2}" 
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(34.1)5fe(p) Re,(p) 

in the Ar spherical or Ar Cartesian coordinates respectively (see also For- 
mulas 25(1.1, 1.2)), for each k — 1, ...,n, 
(34.2) 5"^(p) 5"^ := ^r...^;^", 

(35) As :=E|i|<aa,'S^(p)- 
Then we have the integral formula: 

(36) A5^"(/xQ»,«;p;C) = V/(t)[A5exp(-«(p,t;C))]rft 

in accordance with 1(7) and 2(4). Due to §28.3 the operator A5 has the 
inverse operator for A„+i almost all (po, ■■■■,Pn) in R""^^. Practically, its cal- 
culation may be cumbersome, but finding for an integral inversion formula 
its kernel is sufficient. In view of the inversion Theorem 6 or Corollary 6.1 
and §§19 and 20 we have 

(37) (27r)-" /r„ exp(-M(a + p, t; Q) exp{u{a + p, r; C))dpi...dpn = S{t; r), 
where 

(38) [5,f){T)=J^„f{t)5{t;T)dh...dtr,^f{T) 

at each point r e R", where the original /(r) satisfies Holder's condition. 
That is, the functional 6{t; r) is Ar linear. Thus the inversion of Equation 
(36) is: 

(39) / ( / f{t)xQn{t){[As exp(-u(p+a, t; (Mip+a, t, r; Q}dt)dp^...dpn = /(r), 
SO that 

(40) [A5exp(-M(p+a,t;C))]C(p+a,^,T;C) = (27r)-" exp(-M(p+a, t; C)) exp(-M(p+ 
a,r;C)), 

where the coefficients of A5 commute with generators ij of the Cayley- 
Dickson algebra Ar for each j. Consider at first the alternative case, i.e. 
over the Caylcy-Dickson algebra A,- with r < 3. 

Let by our definition the adjoint operator AJ be defined by the formula 

(41) A*gr]{pA-X) = J2\j\<aSi*S^r]*{p,t;() for any function : A x R" x 
Ar ~> Ar, where t e R", p and C G A; S^r]*{p,t; Q ■= [S^viP,t;C)]*- Any 
Caylcy-Dickson number z G A^j can be written with the help of the iterated 
exponent (see §3) in Av spherical coordinates as 

(42) z = |z|exp(-M(0,0;^)), 

where v > r, ip E Av, u G Ay, Re{ip) = 0. Certainly the phase shift operator 
is isometrical: 

(43) \T^\..T^^;^z\ = \z\ 

for any ki,...,kn G R, since \exp{~u{0,0; Im{tp))\ = 1 for each ijj G Av, 
while T^...T^"e-"(°'°;-f"*W) - exp{-u{0,0; Im{ij) - {kih + ... + /c„«„)7r/2)} 
(see §12). 
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In the Ar Cartesian coordinates each Cayley-Dickson number can be pre- 
sented as: 

(42.1) z = \z\ exp((/)M), where G R is a real parameter, M is a purely 
imaginary Cayley-Dickson number (see also §3 in pTj [T6]). Therefore, we 
deduce that 

(44) \Asexp{-u{p+a,t;C))\ = exp{-{po+a)si-Co)\As exp{-u{Im{p),t; Im{C)))\, 
since R is the center of the Cayley-Dickson algebra Av and po, a, Co, si G R, 

Si = Si(t), where particularly A5I := A5e~"^'''°'^'*|(j=o (see also Formulas 
12(3.1-3.7)). 

Then expressing ^ from (40) and using Formulas (41, 42, 42.1, 44) we infer, 
that 

(45) aP,t,T;0 = (27r)-[AJexp(-M(Jm(p),t; Jm(C))] 
[exp{-u{Im{p),t; Im{C)) exp{u{p, r; C))]\Asexp{-u{Im{p),t; /m(C))|"^ 

since = z* j\z\^ for each non-zero Cayley-Dickson number z G f > 1, 
where Im{p) = piii + ...+ pJu, P = Wo + ••• + Pnin, Po = Re{p). 

Generally, for r > 4, Formula (45) gives the integral kernel C,{p,t,T]() 
for any restriction of ^ on the octonion subalgebra alg-R^Ni, N2, N4) embed- 
ded into Ar- In view of §28.3 C, is unique and is defined by (45) on each 
subalgebra alg^^Ni, N2, N4), consequently. Formula (45) expresses ^ by all 
variables p,^ E Ar and t, t G R". Applying Formulas (39,45) and 28.2(A) 
to Equation (34), when Condition 8(3) is satisfied, we deduce, that 

(46) (/XQ")(r) = / (/ g{t)xQ^it)[exp{-u{p+a,t;C))^ip+a,t,T;C)]dt)dpi...dpn- 

JR" JR" 

E a, E (-1)' 

lil<" l<|(/j)|, 0<mfc, 0<qk, hk=sign{lkjk); mk+qk+hk=jk; 9fc=0 for lkjk=0, V fe=l,...,n; {0e{0,l,2}" 

/ (/ [d^'^^f{t^''ydtf..Ml-][{S"'{p)exp{-uip+a,t^'^^-Xy)^^^^ 

where dimjidQ'^i-^ = n — \h{lj)\, t^^^^ G dQ'li-) in accordance with §28.1, 5"'(p) 
is given by Formulas (34.1,34.2) above. 

For simplicity the zero phase parameter ^ = in (46) can be taken. In 
the particular case = R" all terms with Jqq^ vanish. 

Terms of the form /Rn[{5™(p) exp(— M(p-|-a, t; C))}^{p+C'^ ^5 0]dpi-..dpn 
in Formula (46) can be interpreted as left Ar linear functionals due to Fubini's 
theorem and §§19 and 20, where 5° = /. 

For the second order operator from (14) one gets: 

(47) As = (ELi ^h[Sh{p)?) + PnSnip) + 00 and 

(48) {fxu)it)=f if g{t)xu{t)[exp{-u{p+a,t;Cm{P,t,r;C)]dt)dpi...dpn- 

JR" JR" 

/ (/ /(t')[{(/9(0 + P(^',P))exp(-w(p+a,t;C))K(p,t',r;C)]cit')c?Pi-c?Pn- 

JR" JdUo 
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/ ( / a{t'){df{t')/dv)[BM-u{p + a, t- C))e(p, t\ T- C)]dt')dp^...dprr. 
J-R" JdUo 

For a calculation of the appearing integrals the generalized Jordan lemma 
(see §§23 and 24 in [18]) and residues of functions at poles corresponding to 
zeros I A5 exp(— M(/m(p), t; Im{())) \ = by variables pi, ...,Pn can be used. 
Take g{t) = 6{y; t), where y G R" is a parameter, then 

(49) / (/ 6{y;t)[exp{-u{p + a,t;())C{p + a,t,T;C)]dt)dpi...dpn 

JR" JR" 

= / bwi-u{p + a,y;C))^iP + a,y,T;C)]dpi...dpn=: £{y;T) 

JR" 

is the fundamental solution in the class of generalized functions, where 

(50) AS{y;t) = 6{y;t), 

(51) I^n6{y;t)f{t)dt = f{y) 

for each continuous function f{t) from the space L^(R'^, Ar)] At is the partial 
differential operator as above acting by the variables t = {ti, ...,tn) (see also 
§§19, 20 and 33-35). 

29. The decomposition theorem of partial differential operators 
over the Cay ley- Dickson algebras. 

We consider a partial differential operator of order u: 

(1) Afix) = J2 a„(x)9°/(x), 

\a\<u 

where 9"/ = d^"'^ f{x) /Oxq" ...dx'^" , x = xoio + ...Xnin, xj E R for each j, 
l<n = 2^ — 1, a = (ao, a„), |a| = ao + ... + an, < aj E Z. By the 
definition this means that the principal symbol 

(2) Ao := E aa(x)5° 

\a\=u 

has a so that |a| = u and 3ia{x) G Ar is not identically zero on a domain 
U in Ar- As usually C''{U,Ar) denotes the space of k times continuously 
differentiable functions by all real variables xq, ...,Xn on U with values in Ar, 
while the x-differentiability corresponds to the super-differentiability by the 
Cayley-Dickson variable x. 

Speaking about locally constant or locally differentiable coefficients we 
shall undermine that a domain U is the union of subdomains satisfying 
conditions 28{Dl,i - vii) and n U^' = dW n dU'' for each j ^ k. All 
coefficients a^ are either constant or differentiable of the same class on each 
IntiU^) with the continuous extensions on . More generally it is up to a 
C" or x-differentiable diffeomorphism of U respectively. 

If an operator A is of the odd order u = 2s — 1, then an operator E of 
the even order m + 1 = 2s by variables {t, x) exists so that 
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(3) Eg{t,x)\t=o = Ag{0,x) for any g e C"+^([c,d] x U,Ar), where t E 
[c, d] C R, c < < d, for example, x) = d{tAg{t, x))/dt. 

Therefore, it remains the case of the operator A of the even order u — 2s. 
Take z — zoio + ... + -22" -1^2"-! £ Av, zj e R. Operators depending on a less 
set zi^, zi^ of variables can be considered as restrictions of operators by all 
variables on spaces of functions constant by variables Zs with s ^ {^i, In}- 

Theorem. Let A — A^ be a partial differential operator of an even 
order u — 2s with locally constant or variable or x-differentiable on U 
coefficients 3ia{x) € Ar such that it has the form 

(4) Af = c„,i(S„,i/) + ... + Cu^k{Bu,kf), 'where each 

(5) -Bu,p = Bu^pfi + Qu-i,p 

is a partial differential operator by variables Xm^ ,^+,„+rn„^^_^+i,...,Xm^ ,^+...+rnu,p 
and of the order u, rriufi — 0, Cu^k{x) e Ar for each k, its principal part 

is elliptic with real coefficients SLp^2a{x) ^ 0, either < r < 3 and f e 
C^{U, Ar), orr>A and f e C^{U, R). Then three partial differential opera- 
tors and Tf and Q of orders s andp withp <u — l with locally constant or 
variable or x-differentiable correspondingly on U coefficients with values 
in Av exist, r < v, such that 
(7) Af^T^{Tlf) + Qf. 

Proof. Certainly we have ordQu-i,p < u — 1, ord{A — Aq) < u — 1. We 
choose the following operators: 

(8) T7(x) = E E (9"/(x))KV'p,a] and 

l"l<s, a(j=OVg<(m„,i+...+mu,p_i+l) and q>{mu,i+---+mu,p) 

(9) T7(a;) = X: Yl (5"/(^))KV^;a], 

^~^\a\<s, aq=OVq<{mu,i+...+mu,p-l+l) and q>{mu.i+---+m.n,p) 

where = Cu^p for all p and ippa{x) = —Sip^2a{x) for each p and x, Wp G Ar, 
'ipp,a{x) G Ar,v and ipp^a{x) is purely imaginary for SLp^2a{x) > for all a and x, 
Re{wplm{iljp^a)) = for all p and a, Im{x) = {x — x*)/2, v > r. Here Ar^v = 
Ay/Ar is the real quotient algebra. The algebra Ar^v has the generators ij2r, 
J = 0, 2^-'--!. A natural number v so that 2^-''-l > EJ=i Eg=o (™^+'^"^) 
is sufficient, where = m\/{q\{m — q)\) denotes the binomial coefficient, 
("^"'7^^) number of different solutions of the equation «! + ... + = q 

in non-negative integers aj. We have either 8°"^^ f G for < r < 3 
or (9"+'^/ G R for r > 4. Therefore, we can take ipp^aix) G i2''gR, where 
q = q{p, a) > 1, q{j>^, oi^) 7^ q{;p, a) when {p, a) ^ o^). 

Thus Decomposition (7) is valid due to the following. For h = 8°"^^ f{z) 
and 1 = i2rp and w E Ar one has the identities: 
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(10) {b{wV)){w*l) = {{wb)l){w*l) = -w{wb) = -w% and 

(11) {{{b\)w*)V)w — {{{bw)l)V)w — —{bw)w — —bw"^ in the considered 
here cases, since Ar is alternative for r < 3 while R is the center of the 
Cayley-Dickson algebra (see Formulas (M1,M2) in the introduction). 

This decomposition of the operator is generally up to a partial dif- 
ferential operator of order not greater, than (2s — 1): 

(12) Qf{x) = E|a|<.,|^|<.;.<a,.</3,|.+.|>o[n?:o' fi^) 

where operators and Tf are already written in accordance with the general 
form 

(13) TV(x) = E|a|<.(9"/(x))77a(x); 

(14) Tlf{x)^j:m,{d^f{x))r]}{x). 

When A in (3) is with constant coefficients, then the coefficients Wp and 
t/jp^a for T"* and T^" can also be chosen constant and Q — 0. 

30. Corollary. Let suppositions of Theorem 29 be satisfied. Then a 
change of variables locally constant or variable or x-differentiable on U 
correspondingly exists so that the principal part A2,o of A2 becomes with con- 
stant coefficients, when ap^2a > for each p, a and x. 

31. Corollary. If two operators E — A2S and A — A2S-1 are related by 
Equation 29{?>), and A2S is presented in accordance with Formulas ^5(4,5), 
then three operators , and Q of orders s, s — 1 and 2s — 2 exist so 
that 

(1) A2,_i = T^T^-i + g. 

Proof. It remains to verify that ord{Q) < 2s — 2 in the case of ^25-1, 
where Q — {d{tA2s-i) / dt — T^T\'\\t=Q. Indeed, the form \{E) corresponding 
to E is of degree 2s — 1 by x and each addendum of degree 2s in it is of 
degree not less than 1 by t, consequently, the product of forms A(T5)A(Tf) 
corresponding to and T\ is also of degree 2s — 1 by x and each addendum 
of degree 2s in it is of degree not less than 1 by t. But the principal parts of 
\{E) and A(T5)A(Tf) coincide identically by variables {t, x), hence ord{{E — 
T*T^}|(=o) < 2s — 2. Let a{t,x) and h{t,x) be coefficients from T\ and T*. 
Using the identities 

a{t,x)dtd^tg{x) — a{t,x)d^ g{x) and 

h{t,x)d^dt[a{t,x)dng{x)] = h{t,x)d^[{dta{t,x))dn g{x)\ 
for any functions g'(x) e C^*~^ and a(i, x) e , ord[{h{t,x)d^),{a{t,x)dn)\\i=Q < 
2s - 2, where dt = d/dt, \^\ < s - 1, I7I < s, [A, B] := AB - BA denotes 
the commutator of two operators, we reduce (T^Tf + Qi)\t=o from Formula 
29(7) to the form prescribes by equation (1). 

32. We consider operators of the form: 
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(1) ir^+PIr)f{z) := {Eo<\a\<k{d'^mVa{z)} + fiz)Piz), 

with r]a{z) e Av, a = (ckq, Q;2r_i), < ctfe e N for each k, \a\ — + ... + 
a2r-,, Pirfiz) := f{z)(3, 

d^f{z) := d\^\f{z)/dz^°...dz2rj{\ 2<r<v<oo, P{z) e A, zo, ^2-1 e 

R, 2; = ZqIo + ... + Z2r-li2r-l. 

Proposition. The operator {T'^ + /3)*{T'^ + P) is elliptic on the space 

Proof. We estabhsh the identity 

(2) {ay)z* + {az)y* = a2Re{yz*) 

for any a,y,z & Av It is sufficient to prove Equahty (2) for any three basic 
generators of the Cayley-Dickson algebra Av, since the real field R is its 
center, while the multiplication in Av is distributive {a -\- y)z — az -\- yz and 
{{aa){py)){'yz*) — {aP^){{ay)z*) for all a, ^,7 e R and a,y,z e Av. If 
a — io, then (2) is evident, since yz* + zy* — yz* + (yz*)* — 2Re{yz*). If 
y — io, then {ay)z* -\-{az)y* — az*+az — a2Re{z) — a2Re{yz*). Analogously 
for z — io. 

For three purely imaginary generators ip, is, ik consider the minimal Cayley- 
Dickson algebra $ = o,lg^{ip, ig, ik) over the real field generated by them. If 
it is associative, then it is isomorphic with either the complex field C or the 
quaternion skew field H, so that {ay)z* -\-{az)y* — a{yz* + zy*) — a2Re{yz*). 

If the algebra $ is isomorphic with the octonion algebra, then we use 
Formulas (M1,M2) from the introduction for either a, y e H and 2; = 1 or 
a,z eH and y — I. This gives (2) in all cases, since the algebra alg^{ip, ig) 
with two basic generators ip and is is always associative. Particularly, if 
y — ^ z^ik, s, k>l, then the result in (2) is zero. 

Using (2) we get more generally, that 

(3) i{ay)z*)b* + {{az)y*)b* = {a2Re{yz*))b* = {ab*)2Re{yz*), 
consequently, 

(4) {{ay)z*)b* + {{az)y*)b* + {{by)z*)a* + {{bz)y*)a* = me{ab*)Re{yz*) 
for any Cayley-Dickson numbers a, b,y,z & Av. In view of Formulas (1, 4) the 
form corresponding to the principal symbol of the operator (T'^-|-^)*(T'^-|-^) 
is with real coefficients, of degree 2k and non-negative definite, consequently, 
the operator (T'^ + ^)*{T'^ + (3) is elliptic. 

33. Fundamental solutions. Let either y be a real y = ^ or com- 
plexified Y — {Av)c or quaternionified Y — {Av)u Cayley-Dickson algebra 
(see §28). Consider the space B(R",Y) (see §19) supphed with a topology 
in it is given by the countable family of semi-norms 

(1) p^,k{f) sup,,i,„ 1(1 + \x\)'^d-f{x)\, 
where k — 0,1,2,...; a — (cti, «„), < aj & Z. On this space we 
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take the space B'{Ii"',Y)i of all Y valued continuous generalized functions 
(functionals) of the form 

(2) / = /oio + ... + f2^-\i2^-i and g = qqIq + ... + 92^-112^-1, where 
fj and Qj e B'(R",y), with restrictions on B(R",R) being real or Ci or 
^j,K,L -valued generalized functions fo, f2^-i, go, 92'"-i respectively. Let 
(j) = 00^0 + ... + 02''-i^2''-i with 00, 02"-! ^ B(R"', R), then 

(3) [/, (j)) — Z]fej=o[/j , (f)k)ikij- We define their convolution as 

(4) [/ *g,(f))^ J:f,k~Jo{[fj * 9k, (t>)h)ik for each e B(R", F). As usually 

(5) (/ * g){x) = /(x - y) * = /(y) * g{x - y) 

for all X, y e R" due to (4), since the latter equality (5) is satisfied for each 
pair fj and gk- Thus a solution of the equation 

(6) (T^ + p)f^g in B(R", Y) or in the space B'(R", Y)i is: 

(7) / = Sr^+i3 * g, where Sr^^/s denotes a fundamental solution of the 
equation 

(8) (T^ + p)Sr+^ = 5, (5,0) = 0(0). The fundamental solution of the 
equation 

(9) AoV = 5 with Ao = (T^ + ^)(Tf + ^1) 

using Equalities 32(2 — 4) can be written as the convolution 

(10) V =: Vao = ^r^+is * ^tJi+/3i- 

More generally we can consider the equation 

(11) Af = g with A = ^ + (T2 + (32), 

where Aq = {T + (3)(Ti + f3i), T, Ti, T2 are operators of orders s, si and S2 
respectively given by 32(1) with z-differentiable coefficients. For T2 + /92 = 
this equation was solved above. Suppose now, that the operator T2 + (32 is 
non-zero. 

To solve Equation (11) on a domain U one can write it as the system: 

(12) (Ti + A)/ = (T + P)gi = <7 - (T2 + (32) f. 

Find at first a fundamental solution Va of Equation (11) ioi g = 5. We have: 

(13) / = £^Ti+/3i *gi = ^T2+/32 *ig- + l^)gi), consequently, 

(13.1) ^Ti+/3i * ^1 + ^T2+/32 * (C^ + (3)gi) = ^T2+/32 * 9- 

In accordance with (3 — 5) and 32(1) the identity is satisfied: [6x2+^2 * ((^ + 
(3)gi), 0o) = [(TT + (3){£r2+i32 * 9i), 0o). Thus (13) is equivalent to 

(14) ^Ti+/3i * ^1 + (T + (3){Sr2+l32 * gi) = ^T2+/32 

for g = 5, since Sr2+/32 * ^ = ^T2+/32- 

We consider the Fourier transform F by real variables with the generator 
i commuting with ij for each j = 0, 2" — 1 such that 

(Fl) {Fg){y) = e^'^y'^^ g{x)dx^...dx,, 
for any g G L^(R",^^), i.e. f-^n \g{x)\dxi...dxn < 00, where g : R" — )■ y 
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is an integrable function, {y,x) = Xiyi + ... + XnUn, x = (xi, G R", 

Xj G R for every j. The inverse Fourier transform is: 
(F2) {F-^g)iy) = i27r)~^ e''^y'^^g{x)dxi...dxn. 

For a generalized function / from the space S'(R", Y)i its Fourier trans- 
form is defined by the formula 

(F3) (F/, 0) = (/, (F-7, 0) = (/, F-V). 

In view of (2 — 5) the Fourier transform of (14) gives: 

(15) [F(fx,+ft)][F(^7i)] + E?:oMi^((T + /3)A2+/5j][^(^?i)]^. = ^(^x.+z^J 
for g = 5. With generators io; "^2"-!, "^oi; •••,i2''-ii the latter equation gives 
the linear system of 2^+^ equations over the real field, or 2^+^ equations when 
Y = {Av)ii- From it F{gi) and using the inverse transform a generalized 
function gi can be found, since F{g) = F{go)io + ... + F{g2v^i)i2^-i and 
F~^{g) = F^^{gQ)iQ + ... + F~^{g2v ^1)12^^1 (see also the Fourier transform of 
real and complex generalized functions in [SlEn])- Then 

(16) Va = ^^Ti+/3i * gi and f = Va* g gives the solution of (11), where gi 
was calculated from (15). 

Let vr^ : {Av)-h. (^r)H be the R-linear projection operator defined as 
the sum of projection operators tcq + ... + 7r2r_i, where ttj : {Ay)^ — ?■ Hij, 

(17) Tijih) = hjij, h = hjij, hj G iij^K,L, that gives the corre- 
sponding restrictions when hj G Ci or hj G R for j = 0, ...,2'" — 1. Indeed, 
Formulas 2(5, 6) have the natural extension on {Av)ii, since the generators 
J, K and L commute with ij for each j. 

Finally, the restriction from the domain in A^ onto the initial domain of 
real variables in the real shadow and the extraction of tt^ o / g with the 
help of Formulas 2(5, 6) gives the reduction of a solution from Av to Ar. 

Theorems 29, Proposition 32 and Corollaries 30, 31 together with formu- 
las of this section provide the algorithm for subsequent resolution of partial 
differential equations for s,s — 1, ...,2, because principal parts of operators 
A2 on the final step are with constant coefficients. A residue term Q of 
the first order can be integrated along a path using a non- commutative line 
integration over the Cayley-Dickson algebra fTli [TB] . 

34. Multiparameter transforms of generalized functions. 

If G B(R", Y) and g G i3'(R", Y)i (see §§19 and 33) we put 

(1) E?lo'['^"(^7.;^;p;C),0)^, := E?lo'fe,'^"(0;«;p;C))^, or shortly 

(2) EWl9Je--^''''''^\<|>)^, = E^M^?., ^e-fe*^^)^,. 

If the support supp{g) of g is contained in a domain U, then it is sufficient to 
take a base function with the restriction 0|(7 G B{U, Y) and any 0|Rn\{/ G 

34.1. Remark. It is possible to use Theorem 29, Corollaries 30 and 
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31, Proposition 32 and §33 for solutions of definite differential equations 
with variable coefficients. For this purpose one can present an operator 
A as the composition A = TTi + Q, where ord{A) = ord{T) + ord{Ti), 
ord{Q) < ord{A) — 1, T and Ti are operators with variable coefficients, T*T 
and TJT 1 are elliptic operators with constant coefficients of their principal 
symbols at least. Then use Formulas 33(1 — 16) to find fundamental solutions 
£r, £ti and £a or iterate this procedure (see also §35). A generalization of 
Feynman's formula over the Cayley-Dickson algebras for the second order 
partial differential operators with the first order addendum Q with variable 
coefficients from [20] also can be used. 
35. Examples. Let 

(1) Af{t) = Y.Uid'f{t)/dt])c, 

be the operator with constant coefficients Cj G At, \cj\ = 1, by the variables 
ti,...,tn, n > 2. We suppose that Cj are such that the minimal subalge- 
bra alg-£i{cj,Ck) containing cj and Ck is alternative for each j and k and 
l(-(cfcf)...)cy2| = i. Since 

(2) dfit)/dt, = T.l=i{df{t{s))/ds,){dsu/dt,) = n=idf{t{s))/dsk, the 
operator A takes the form 

(3) Af = j:U{J:i<k,b<j{d'f{tis))/ds,ds,))cj, 

where sj = tj + ... + tn for each j. Therefore, by Theorem 12 and Formulas 
25(50) and 28(6) we get: 

(4) J^"{Af;u;p;0 = E"=i{[Re,(p)]'i^:(p; C)}c, for uip,t;0 either in 
Ar spherical or Ar Cartesian coordinates with the corresponding operators 
Rej(p) (see also Formulas 25(1.1, 1.2)). 

On the other hand, 

(5) J'"(5; u;p; Q = e-"(?''°'^) = e-"(°'°'^) in accordance with Formula 20(2). 
The delta function 6{t) is invariant relative to any invertible linear operator 
C : R"" — >■ R" with the determinant | det(C)| = 1, since 

/ 6{Cx)(j){x)dx= f 6{y)(f){C-^y)\det{C)\dy = (j){C-^0) = (j){0). 

Thus 

(5) -F"(C(A/); u- p- C) = ^"(A/; u; p; Q 
for any Fundamental solution /, where Cg{t) := g{Ct), Af = 6. If C : 
R" — 7- R*^ is an invertible linear operator and ^ = Ct, q = Cp, = CQ, then 
t = C^^C,, p = C~^q and ( = C^^C,' . In the multiparameter noncommuta- 
tive transform J-"" there are the corresponding variables {tj,pj,Q)- This is 
accomphshed in particular for the operator C(ti,...,t„) = (si,...,s„). The 
operator C^^ transforms the right side of Formula (4), when it is written in 
the Ar spherical coordinates, into J2]=i{{po + Qj^e^Y {q] ()}cj. The Cayley- 
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Dickson number q = qo + qiii + ... + g„iri can be written as g = go + qAiM, 
where \M\ = 1, M is a purely imaginary Cayley-Dickson number, qM € R, 
qMM = qiii + ... + g„i„, since go = Re{q). After a suitable automorphism 
6 : Ar ^ Ar we can take ^(g) = go + qMH, so that 6{x) = x for any real 
number. The functions [Sj=i Q'jSg^Cj] and [Z]j=iPjSg^Cj] are even by each 
variable g^ and pj respectively. Therefore, we deduce in accordance with 
(5) and 2(3,4) and Corollary 6.1 with parameters po = and C = and 
Cj G {—1,1} for each j that 

(6) (-F")-Hl/[E"=i{Ei<M<.P'^Se,j9,SeJc,];«;y;C) = -[^?, e^d^l'I'^D) 
in the Ar spherical coordinates, where g = l/[J2]=i or 

(6.1) iJ'T\U[J:U{ppl}c,];u;y;0 = -[^7, e^([^l'W)) 
in the Ar Cartesian coordinates, where g = l/[Ej=i PjCj], N = y/\y\ for 
y ^ 0, N = ii hi y = 0, y = yih + ... + G [y] = ivu Vn) e R", 
(M, [q]) = E"=i yjqj, since S^^ cos(0 + Cfc) = cos(0 + + tt) = - cos(0 + Cfc) 
and Sg^ sin(0 + (k) = sin(0 + Cfc + tt) = — sin(0 + Cfc) for each k. 

Particularly, we take cj = 1 for each j = 1, fc+ and cj = —1 for any 
j = fc+ + 1, ...,n, where 1 < A;_|_ < n. Thus the inverse Laplace transform for 
go = and C = in accordance with Formulas 2(1 — 4) reduces to 

(7) (-F")-Hl/[E"=i{Ei<M<,PfcSe,P6SgJc,];M;y;C) = 

(27r)-" exp(i(gi|/i + ... + g.y„))(l/[E -ii g| - E^=fc^+i g^])f^gi...rfgn 
in the .4^ spherical coordinates and 

(7.1) ij^-)-\i/[j:upplM^^-^y^O = 

(27r)-"/i,„exp(i(pi2/i + ... +p„l/„))(l/[E-;iPf - E-=fc^+iP|])t^Pi-t^Pn 
in the Ar Cartesian coordinates, 

since for any even function its cosine Fourier transform coincides with the 
Fourier transform. 

The inverse Fourier transform {F~^g){x) = {2TT)~'^{Fg){—x) =: \E'„ of the 
functions g = l/{E]=iz]) forn > 3 and P(l/(E?=i4)) forn = 2 in the 
class of the generalized functions is known (see [5] and §§9.7 and 11.8 [29] ) 
and gives 

(8) ^„(zi,...,z„) = C„(E^=i4)'-"/' 3 < n, where C„ = -l/[{n - 
2)cr„], a-a = 47r"/^/r((r2/2) — 1) denotes the surface of the unit sphere in R", 
r(x) denotes Euler's gamma-function, while 

(9) ^2{zi, Z2) = C2 ln(E?=i z]) for n = 2, where C2 = l/(47r). 

Thus the technique of §2 over the Cayley-Dickson algebra has permitted to 
get the solution of the Laplace operator. 
For the function 

(10) P{x) = Ejii ~ Ej=fc++i ^'j with 1 < k+ < n the generalized func- 
tions {P{x) + iO)^ and {P{x) — iO)^ are defined for any A G C = R © iR (see 
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Chapter 3 in [S]). The function has the cone surface P{zi, = of 

zeros, so that for the correct definition of generahzed functions corresponding 
to the generahzed functions 

(11) {P{x) + ciO)^ = hmo<c.,e^o(^(a:)^ + e^f/"^ Qx.^{\\arg{P{x) + ice)) 
with either c=— lorc=l were introduced. Therefore, the identity 

(12) F(vl/,^,„_,J(x) = -{Y!;Ux] - E-=fc^+ix2)[F(v^,^,„_,J(x)]2 or 

(13) F{m) = -l/(P(x) + ciO) follows, where c = -1 or c = 1. 

The inverse Fourier transform in the class of the generalized functions is: 

(14) F^\{P{x) + ciOy){zi, Zn) = exp(-7rc(n-A;+)i/2)22^+"7r"/2r(A + 
n/2)(g(zi,...,^0-ciO)-^-"/2)/[r(-A)|D|V2] 

for each A G C and n > 3 (see §1V.2.6 [5J), where D = det{gj^k) denotes 
a discriminant of the quadratic form P(^x) = X/?a;=i 9j,kXj-^ky 

while Q{y) = 

Yl]',k=i 9^'^XjXk is the dual quadratic form so that I]fc=i 9^'''9k,i = for all 
j, I; 6j = 1 for j = I and Sj = for j ^ I. In the particular case of n = 2 the 
inverse Fourier transform is given by the formula: 

(15) F-i((P(x)+ciO)-i)(zi,22) = -4-i|D|-V2exp(-7rc(n-fc+)i/2)ln(Q(zi, 
ciO). Making the inverse Fourier transform of the function —1/ (P(x)+iO) 
in this particular case of A = —1 we get two complex conjugated fundamental 
solutions 

(16) '^k+,n-k+izi, ...,Zn) = -exp{7ic{n-k+)i/2)T{{n/2)-l){Q{zi, ...,z„) + 
ci0)i~("/2)/(4^n/2) for 3 < n and 1 < k+ < n, while 

(17) ^i,i(;zi,;z2) = 4^^exp(7rc(n - k+)i/2)\n{Q{zi, Z2) + ciO) for n = 2, 
where either c=lorc=— 1. 

Generally for the operator A given by Formula (1) we get P{x) = Po{x) + 
Pi{x), where Po{x) = Yl]=ix'jRe{cj) and Pi{x) = YIj=ix']Im{cj) are the real 
and imaginary parts of P, Im{z) = z — Re{z) for any Cayley-Dickson number 
z. Take 1 = and consider the form P(x) + eel with e 7^ and either c = 1 
or c = —1, then P(x) + eel 7^ for each x G R". We put 

(18) (P(x) + doy = \imo<ce,e^o{P{xf + e^y/^ exp{iXArg{P{x) + Ice)). 
Consider A G R, the generalized function (P(x)^ + e^)'*'/^ exp(iAylr(7(P(x) + 
Ice)) is non-degenerate and for it the Fourier transform is defined. The limit 
limo<c€,e^o gives by our definition the Fourier transform of {P{x) + clO)^. 
Since 

(19) Cj{(3j + J2l<k<n,kytj Cj '^Ckl^k) = Yl^=l 

for all Pj G R and any 1 < j < ri in accordance with the conditions imposed 
on Cj at the beginning of this section and iNj = Nji for each j, the Fourier 
transform with the generator i can be accomplished subsequently by each 
variable using Identity (19). The transform Xj 1— )■ (cjY^'^Xj is diagonal and 
|(...((c|;^^C2^^)...)cy^| = 1, so we can put \D\ = 1. 
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Each Cayley-Dickson number can be presented in the polar form z = 
\z\e'^^^ , (j) G R, 101 < TT, Misa purely imaginary Cayley-Dickson num- 
ber \M\ = 1, Arg{z) = {</) + 27ck)M has the countable number of val- 
ues, k E 2i (see §3 in [TTl US]). Therefore, we choose the branch z^^"^ = 
\z\^/'^exp{{Argz)/2), \z\^/'^ > for ;z 7^ 0, with \Arg{z)\ < tt, Arg{M) = 
Mtt /2 for each purely imaginary M with \M\ = 1. 

We treat the iterated integral as in §6, i.e. with the same order of brackets. 
Taking initially cj G R and considering the complex analytic extension of 
formulas given above in each complex plane R © A'jR by Cj for each j by 
induction from 1 to n, when Cj is not real in the operator A, Im{cj) G UNj, 
we get the fundamental solutions for A with the form (P(a;) + clO)'*' instead of 
{P{x)+ciO)^ with multipliers (...(ci^^C2^^)...)c^/^ instead of exp(7rc(ri — A;+)i/2) 
as above and putting \D\ = 1. Thus 

(20) xl>(^i,...,^0 = -r((n/2)-l)(P*(^i,...,^„)-dO)Mn/2)[(...(cfcf )...)c^^^^ 
for 3 < n, while 

(21) ^(^1, Z2) = ^-\cTcTX'Ln{P\z^, Z2) - clO) for n = 2, 
since c* = cj"^ for \cj\ = 1, yjqj = yjic'I'^'YqjC-'^ ^ while 

{...{dci'^qxdcfq2)..:)dci'^q^ = dqi...dqn[{...{c^^cj'^)...)c'il'^] and 
l(-(cfcf)...)</2| = l. 

36. Partial differential equations with polynomial real coeffi- 
cients. Let 

(1) A = E|a|<m«a(g)5°, 

(^a{(l) = J2i3 o.a,i3(l^ , '■= Qi^-'-Qn", <^a,i3 and / have values as in §28, and 
Af be an original. Using the transform in the Ar Cartesian coordinates we 
take qj = tj for each j, while using the transform in Ar spherical coordinates 
we choose qj = Sj (t) for each j . Then 

(2) T-iAf;u;p;C) = Ef^i-iy^%ip)d^p 

a^AlPo + PiSeJ^ip? S-^..p^S^„")]F"(p; C) = G^ip; () 
in the Ar spherical coordinates and 

(2.1) T-iAf; u;p; Q = Epi-iy^%ip)d^p 

(E/3a«,/3[pO+PlSeJ°4PO+P2Se,]"^..[Po+PnSeJ"")F-(p; C) = ^"(p; C) 

in the Ar Cartesian coordinates (see Theorems 12 and 13 above). It may 
happen that the second differential equation is simpler than the initial one: 

(3) Af = g. 

For example, when coefficients depend only on one variable then the 
second differential equation is ordinary and linear. 

37. Noncommutative transforms of products and convolutions 
of functions in the Ar spherical coordinates. 

For any Cayley-Dickson number z = zqIq + ... + Z2r^ii2r-i we consider 
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projections 

(1) ej(z) = Zj, Zj e R or Ci or lij^K,L, j = 0, .., 2^ - 1, Ojiz) = 7r,(z)i*, 
given by Formulas 2(5, 6) and 33(17). We define the following operators 

(2) 7e.,,-(F"(p; 0) := F"(po, {-iT^Pi, (-l)'^^+-^^-p,+i-5,.„,P,+2-5,,„, 

-,Pn; Co, (-l)"iCl+™l/2, (-l)"^+^-*-"0+l-5,,„+7rQ;,-+i_5.„/2, 0+2-5,,„, Cn) 

on images F", 2^"^ < n < 2^ - 1, j = 0, For and pj e {0, 1} their 
sum aj + Pj is considered by {mod 2), i.e. in the ring Z2 = Z/(2Z), for two 
vectors a and ^ e {0, 1}^"^"^ their sum is considered componentwise in Z2. 
Let 

n V-l 

(3) J-"(/;i.;p;C) = E E ^.■(•F"(^fe(/); OK^i, 

also Fj^(p; C) := Efc=o^ 0)4 for an original /, where u{p, t; Q 

is given by Formulas 2(1,2,2.1). If / is real or Ci or tlj^K,L -valued, then 

Theorem. If f and g are two originals, then 

(4) T^ifg;u;p;0 = E -=0 Ea,/36{o,i}«(-l)"^+^^'-'^+^-H^a,,(i^nP-?o; C- 
V)) * {T^/3,j{G]{p + qo- po; r]))ij, 

(4.1) J^-{f*g;u;p;0 = E^=o Ea,^e{o,i}"(-l)"^+^^'"'^+^-n^aj(^"(p; C - 
r])){n^AG]{p;v))ij, 

whenever J^"'{fg), J^"(/); J^'"{g) exist, where 1 <n < 2^ — 1, 2 < r; ak+/3k — 
1 {mod 2) for k < j or k — j + 1 — n, ak + (3^ — {mod 2) for k — j + 1 < n 
and ak — — for k > j + 1 in the j-th addendum on the right of Formulas 
(4,4.1); the convolution is by (pi, ...,Pn) in (4), at the same time qo and 
T) & Ar are fixed. 

Proof. The product of two originals can be written in the form: 

(5) f{t)g{t) = E?Lo^ E.,. ek{f{m{g{t))ij. 

The functions 9k{f) and 9i{g) are real or Ci or 11j^k,l valued respectively. 
The non-commutative transform of fg is: 

(6) -^"(/^)(p;C)= / f{t)g{t)eM-<p,t;C))dt^ 

{ / {f{t)g{t))e-P°'' cos(piSi + Cl)^oc^^}+ 

n— 1 „ 

{E / sin(piSi+Ci)... sin(pj_iSj_i+Ci-i) cos{pjSj+Q)ij-idt}+ 

/ (/(t)5(t))e-f°^i sin(pisi + Ci)... sin(p„s„ + Qn)indt. 

7R" 

On the other hand, 

(7) / /(t)^(t)e-^°'^+'^'=>^'^+^^'^^^di = 

^R" 
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Jr." jr" jr" 

where k = 1,2, ...,n, G { — 1,1}. Therefore, using Euler's formula e"^ = 
cos(0) + isin(0) and the trigonometric formulas cos(0 + ■?/') = cos(0) cos(-(/^) — 
sin(0) sin(-?/'), sin(0 + -i/;) = sin(0) cos('?/') + cos(0) sin(-?/') for all (f), ip & 
R, and Formulas (6,7), we deduce expressions for 9j{J^'^{f g)). We get the 
integration by gi, ...,qn, which gives convolutions by the pi, ...,Pn variables. 
Here go £ R- and t] E Ar are any marked numbers. Thus from Formulas 
(5 - 7) and 2(1,2,2.1,4) we deduce Formula (4). 
Moreover, one certainly has 

(8) / {f * g){t)e'''°''^'^U^'''''+^'^^' dt ^ 

JR" 

JR" JR" 

for each 1 < k < n, 7^- G { — 1,1}, since Sj{t) = Sj{t — r) + Sj(r) for all 
j = l,...,n and t, t E R". Thus from Relations (6,8) and 2(1,2,2.1,4) and 
Euler's formula one deduces expressions for 6j{J^"'{f * g)) and Formula (4.1). 

38. Moving boundary problem. 

Let us consider a boundary problem 

(1) Af = g in the half-space t„ > (/)(tn), where 0(0) = and 0(t„) < t„ 
for each < G R. Suppose that the function t„ — 0(t„) =: ip(tn) is 
differentiable and bijective. For example, if < w < 1 and 0(t„) = vtn, 
then the boundary is moving with the speed v. Make the change of variables 
Vn = ^(^n), yi = ti,. ..,?/„„! = t„_i, then tn = ^"Hl/n) and dtn = dSn = 
{dtn/dyn)dyn and due to Theorem 25 we infer that 

(2) r^i Y: h^dfxy.>om;p; C) = E ^A,a. - 1) 

|a|<m |a|<m,0<gn<Q:n — 1 

(Po+SeiPi)>r ■■■P^l^p""-'""'S„,,_(,„+i)e„ J--^'^"(9t«;(y), ^(^ (?/"); C);p; C) 

+ E ha{po + Se,PirpT---PTSa-a,e,J'^Xyn>o{yMyy,p-,0 = G''{p-,C) 
\a\<m 

in the Ar spherical coordinates and 

(2.1) J-( E aa9rx.„>o/(t);p;C)= E a«('^o,a„-l) 

|a|<m \a\<m,0<qn<ctn — i 

{PO + Se,Pir{po+P2Se,r...{PO+Pn-lSe„.,r''-'{PO+PnSe„r--'--' 

J^-^>r{df:w{y),u{p, (i/");C);p;C) 

+ E ^a{P0+Se,Pir{P0+P2Se,r.4P0+PnSeX-J'''{Xyr.>0{y)w{y);P;C) = G^p-Q 

\a\<m 

in the Ar Cartesian coordinates, where w{y) := f{t{y)){dtn/dyn). 
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Expressing J^'^{Xy„>o{y)w{y);p;C) through G'^{p;() and the boundary 
terms J^~^''^" {df^w{y),u{p, (y"');C);P;C) ^ ill §28.3 and making the inverse 
transform 8(4) or 8.1(1), or using the integral kernel ^ as in §28.5, one gets 
a solution w{y) or f{t) = w{y{t)){dyn{tn) /dtn). 

39. Partial differential equations with discontinuous coefficients. 
Consider a domain U and its subdomains U D Ui D ... D Uk satisfying 

Conditions 28(D1, DA^i—vii) so that coefficients of an operator A (see 28(2)) 
are constant on Int{Uk) and onVi — U \ Int{Ui), V2 — Ui \ Int{U2),...,Vk — 
Uk-i \ Int{Uk) and are allowed to be discontinuous at the common borders 
dVj n dUj for each j — 1, /c. Each function fxuj is an original on [/ or a 
generalized function with the support supp{fxuj) C Uj if / is an original or 
a generahzed function on U. Choose operators with constant coefficients 
on and ^^|/„t(y^.) = 0, where = U, so that A\ui^ = A'^,---, A\if. = 
A^ + ... + A'^,..., A\u — A^ + ... + A''. Therefore, in the class of originals or 
generalized functions on U the problem (see 28(1, 2)) can be written as 

(1) Af = g, or 

(2) A^fxv, = gxv„-,A''~^fXVk = gXv^, A'^fxuk = QXu^, 
since xvi + ■■■ + XVk + XUk — Xu- Thus the equivalent problem is: 

(3) A'f = g', A^p = = 

with = fxuk, 9^ = gXUk, also p = fxv^+i, 9^ = 9XVj+i for each 
j — 0,...,k — 1. On dU take the boundary condition in accordance with 
28(5.1). With any boundary conditions in the class of originals or gener- 
alized functions on additional borders dUj \ dU given in accordance with 
28(5.1) a solution p on exists, when the corresponding condition 8(3) is 
satisfied (see Theorems 8 and 28.1). 

Each problem A^ p — g^ can be considered on Uj, since supp{g^) C Uj. 
Extend p by zero on U \Vj for each < j < k — 1. When the right 
side of 28(6) is non-trivial, then p is non-trivial. If p~^ is calculated, then 
the boundary conditions on dU^ \ dU can be chosen in accordance with 
values of p~^ and its corresponding derivatives {d^ f^~^ / dv^)\(QU3\du) for 
some P < ord{A^) in accordance with the operator A^ and the boundary 
conditions 28(5.1) on the boundary dU^ \ dU. Having found for each 
j — 0, k one gets the solution f — p -\- ... -\- on U of Problem (1) with 
the boundary conditions 28(5.1) on dU. 

40. Remark. The multiparameter noncommutative transform over the 
Cayley-Dickson algebras presented above is the natural generahzation of the 
usual complex one-parameter Laplace transform. It opens new opportunities 
for solving partial differential equations of different types. 

It may happen that Theorem 13 is simpler to use, than Theorem 21 



62 



for partial differential equations with real variables. Theorem 13 has an 
advantage that it can be simpler used for partial differential equations of 
complex and hyper-complex variables, because each pair [pi + pjilij) for 
/ 7^ j is the complex variable. In these variants boundary conditions may be 
for F^{p] Q on a, hyperplane Re{p) = a in Ar- 

As it was seen above the appearing integrals are by multidimensional do- 
mains. For their calculations the Fubini's theorem, residues, Jordan Lemma 
and tables of known integrals also can be used. Generally in computational 
mathematics integrals are easier to calculate, than to solve partial differen- 
tial equations numerically. As a rule iterations of algorithms for integrals 
converge faster, than iterations of numerical methods for partial differential 
equations. 

Functions with octonion values may be used to resolve systems of partial 
differential equations. Using conjugations of Cayley-Dickson numbers one 
gets the transition between operators with coefficients either on the left or 
on the right of partial derivatives: [((9"/(x))cq,]* = c* ((9°/(x))*, particularly, 
= for X E R", (9° = d^. 

Using of Formulas 2(5,6) gives variables tj = Zj for z E Ar- So one can 
consider a class of super-differentiable originals f{z),zEVcAr- In the 
class of piecewise on open subsets super-differentiable originals f{z),zE V C 
Ar, with tj = Zj for each j = 1, n, n = 2*" — 1, in the fixed z- representations 
we get the noncommutative transform for f{z)xv{^) relative to the Cayley- 
Dickson variable z E Ar- Therefore, the results given above transfer on this 
variant also. 

Theorem 17 also opens new opportunities to investigate and solve certain 
types of nonlinear partial differential equations using previous results on 
spectral theory of functions of operators [22] • For example, analytic 
functions q{z) in Theorem 17 permit to consider nonlinear operators qic), 
where crf{z) := YlJ=o{df{z)/dzj)ij- It is planned to study in the next paper. 

Partial differential equations with periodic g and / with vector period 
corresponding to may be considered also. Certainly others classes of 
smoothness, for example, Sobolev's or generalized functions can also be con- 
sidered. It is planned in a next paper to consider this and also problems 
with boundary conditions as well as with non-constant coefficients in more 
details. 

The technique described above permits to consider partial differential 
equations of different types and write their solutions in integral forms. If 
appearing integrals can be calculated in elementary or special of generalized 
functions, then this gives the explicit formulas in terms of known functions. 
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In conjunction with the hne integration over the Cayley-Dickson algebras it 
permits to solve some types of non linear partial differential equations. The 
multiparameter Laplace transform over the Cayley-Dickson algebras takes 
into account the boundary conditions. It naturally means the treatment of 
systems of partial differential equations due to the multidimensionality of the 
Cayley-Dickson algebras. 
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